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CHAPTER - ONF 


INTRODUCTION 

Operational research is the application of the methods of science to complex 
problems arising in the direction and management of large system of men, machines, 
materials and money in industry business, government and defence. The distinctive 
approach is to develop a scientific model of the system in corporation measurement of the 
factors such as chance and risk. With which to predict and compare the outcome, 
alternating decisions, strategies, control. The purpose is to help management to determine 
its policy and action scientifically. Operational Research is concerned with scientifically 
deciding how to best design men machine system. Usually under condition requiring the 
allocation of scarce resources. Operational research can be extremely useful in solving 
problems in developing countries as well as developed countries. The developed countries 
if confronted with difficult problem can afford in some money in projects money and 
resources in projects of obscure value, but this is not so far the developing courtiers. As 
they have to be extremely careful in utilizing their resources for getting the optimum results 
which needs the technique of operational Research. The essence of operational research 
activity lies is the construction and use of model. 

A model is a simplified representation of something result. There are many 
conceivable reason why one might prefer to deal with a substitute for the real thing rather 
than with the thing itself. Often motivation is to save the money, time and some other 
valuable commodity. Some time it is to wide risk associate it with tempering of with real 
objects. Sometimes the real environment is so complicated that a represantive model is 



needed just to understand and to communicate with others about it. Queueing analysis has a 
prominent place among various branches of operational research as it deals with 
phenomena that we all encounter as an everyday lives. In fact one is almost encounter some 
form of queue or waiting line in every waking hour. For example such as standing in line in 
ticketing window. There are many subtle cases of queuing behaviors that we experience 
further while making a telephone call. There is a brief delay as a free communication 
channel is to be located while driving a car, a traffic congestion encountered on the open 
highway of the Queueing behavious. The flow of the material through a manufacturing in 
various processing stages represent other form of queue. Queues are found thr oughout 
business in administrative office, medical science, hospital administration, chemical 
industries and in communication technologies. Consider the following example of queues 
of a food processing company may have truck loaded with perishables lined up for over 24 
hours waiting to be unloaded not only spoilages quite only possible perishables already 
loaded. But also those for fields waiting to be loaded of on returning trucks. Another queue 
is found in business frequently associative with inventory. In administration office the 
paper flows throughout the offices and it will be recognized that queue can block this 
communication system. Queue also exist in physical being specific genetic defects in bom 
errors in metabolism in one such disorder a missing principle enzyme cause a substance to 
accumulate wait in queue in blood stream. This eventually damages the brain. Scientist in 
studying the cause and cure of such diseases as cancer and leukemia have found that the 
breaking down of all these process leads to the queue. That makes the disease system. 

In today world it is almost impossible to point out the area in which queue does not 
exist because in the advancement in science and technology has increased the dimension of 
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Ac man and society as a whole to great extent where as the resource are finite and which 
are to be taken into account and therefore there are congestion situation. Queueing theory 
found almost all area demonstrating the divergent nature of the problem. It is applicable to 
inventories best office backlog, large backlog in our court system, catalyst in chemical 
reaction, filtration process, semiconductor noise hospital and the demand for medical care 
nervous reaction such as psychological stress and strain neurosis production line, inspection 
Hne conveyor built. Then there are those barber shop if all short utilities mass 
transportations and communication and so many other facilities. So the life present a very 
large number of queues some are plants or some unexplored consider the amount of time 
wasted while waiting one week, one month or one year period. Hence queues are important 
to managed properly because it is a problem us all directly or indirectly when designing a 
new system or revamping the old one engineers and businessman must be cognizant of 
inherent line. Super market must be planned for adequate check out system and servicing 
system as dissatisfaction may restrict buying or loss of customers. Hospital administrators 
and architects often engage in facilities planning for facilities. Waiting line is common 
place in hospital and they create emotional hardships for the patients waiting for foods, 
admittance, surgery and laboratory work etc. This creates overall frustration, which is to be 
tackled sympathetically. In chemical industry which work on constantly with filtration 
problems must provide that the particles that are gathered by the filter can be cleared at a 
rate greater than the partial arrival rate otherwise the filter will clog and cause the whole 
system to be idled. 

Queueing theory is a mathematical study of waiting line. A basic Queueing model 
has two distinct part. 
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An incorrect design or operating philosophy will magnify all problem areas associated 
with waiting. Thus the study of Queueing theory must incorporate both physical and cost 
characteristics associated with various system. Queues are frequently encountered in 
various networks such as telecommunication, production jobs shops, flexible manufacturing 
process, computer system, urban services system, communication networks, health care 
system, and transportation system. 

A network is a system of lines or channels connecting different points. Some 
example of network are communication line, rail board network, pipeline system, shipping 
line and aviation networks and all those networks some specified commodities is sent from 
certain supply point to certain one because all this networks having multiple resources to be 
taken into the account. For example computer network deal with C.P.U., channels 
memories, communication circuits etc. A queues unavoidable at these facilities these 
complex structure leads to a Queueing networks rather than a single queue with single 
server. Today the information processing industry is one of the fastest growing, most 
dynamic and most glamoures industry on the scientific seen. It highly depends on the 
enormously grown complex networks. The structure and sophistication of these networks 
varies over a considerable range from highly specialized networks design to handle to 
specific task. In a carefully environment to more generalized networks that handle a variety 
of task in highly structured environment. 

Queueing networks are also classified according to service time distribution either 
being exponential or general service time. All customers may belong to the same class with 
the same switching process and requesting similar services or there may be different types 
of costumers. Many practical systems can be modeled of queues open queuing networks 



accepts and losses costumers from to the out side world. Thus the total number of costumer 
in open networks varies with time. 

TANDEM -QUEUES 

Queueing system with several service facilitation in service is known as tandem Queueing 
system. This type of Queueing networks can represent many process of interest in 
manufacturing system, computer system tale communication etc. The co mmuni cation 
system can be typically modeled as networks of interconnected queues by viewing the 
massage as customer communication buffer as waiting line and all activities in necessary 
transmission of the message as services. 

This sort of synchronization has advantages of conceptual simplicity and great generality. 
This leads communication network to view as tandem Queueing system. In the tandem 
Queueing networks there might be an infinite or finite queue capacity at each service unit of 
the network. The tandem queue can be representative for modeling a practical multi stage 
service system in manufacturing process and communication networks. 

It is important to study these queues, for the queues in these areas effect the 
efficiency of the system concerned and cause a huge wastage of precious resources, which 
is ought to be reduced by optimum utilization of resources. 

The production system consists of a number of service facilities where item receive 
service. These system can be very simple or can be very complex with several machines 
linked together in an arbitrary manner, with unequal processing times, limited space for 
temporary storage of items and general routing patterns for different jobs going through the 


5 


system. This is especially true in flexible manufacturing system. These system can also be 
represented as tandem queues by representing the machines as service nodes and the item as 
customer. 
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1.1 REVIEW OF LITERATURE 

The modelling of computer system and data transmission system opened the way, in 
the sixties, to studies of queues characterized by complex service disciplines and have 
created the need of analyze inter connected systems progress in this area has been rapid and 
industrial application have been widely accepted since the seventies at present in the 
computer industry, Queueing network models have resulted in software packages for the 
automatic solution of problems arising in the design of new computer systems in the 
evaluation and improvement of existing system. 

Tandem queues can be representative for modelling practical multistage service system in 
manufacturing process and communicational systems, which attracted many research 
workers to study tandem queue and a lot of work has been done on this line to present a 
brief historic account fallowing works are worth mentioning. 

O.Brien, G.G. (1954) studied the case of two channels in series in the steady state 
with Poisson input and exponential service times and gave expected number and the 
expected waiting times. Burke, P. J. et.al (1956) found the output of a steady state queue, 
independent of the type of service discipline used. Jackson, J.R. (1957), both permitted to 
Poisson arrival to each phase, both from the system and from out side. The service 
distributions are exponential and each phase consists of several parallel channels. An 
important questions to determine in studying tandem queues in the distribution of out put 
from one channel, which the comparises the input into the subsequent channel Reich (1957) 
discussed this problem. Finch, P.D. (1959) shown that it is only when an infinite size queue 
is tolerated that the Poisson result for the output is current. He also gave that N= oo and 
exponential service as a necessary and sufficient condition for the independence of the 
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interdeparture intervals and the independence of the queue length left by a departing unit 
from the interval since the previous departure. 

Friedman (1965) worked on reduction methods for tandem queues. Hillier and 
Boling (1967) studied the tandem (or serial) networks with exponential servers, the most 
basic structural configuration. Neuts (1968) considered two queues in series with a finite 
intermediate waiting room. Fraker (1971) gave approximate techniques for the analysis of 
tandem Queueing system. Tembe et.al (1974) provided the optimal order of service in 
tandem queues. Kobayashi (1974) applied diffusion appromimations to Queueing networks. 
Chandy et.al (1974) did approximate analysis of general Queueing network. Rubin (1974) 
provided path delays in communication networks. Baskett et.al (1975) studied open, closed 
and mixed networks of queues with different classes of customers, Lavenberg (1975) 
discussed stability and maximum departure rate of certain open queuing networks having 
finite capacity constraints. Gelenbe et.al (1976) studied the behaviour of a single queue in a 
general Queueing network Schweitzer (1976) obtained maximum through put in finite 
capacity open queuing network. Kleinrock (1878) studied the communication networks as 
tandem queue. They considered the independent assumption among the service and arrival 
processes. Labetoulle et.al (1977) modelled packed switching communication networks 
with finite buffer size at each node. Dattatreya (1978) shown that the through put of a given 
system is the same as that of the reversed system, and studied tandem queuing system with 
blocking. Harrison (1978) applied diffusion approximation for tandem queues in heavy 
traffic. 

Chandy et.al (1978) gave approximate methods for analyzing Queueing network of 
computing systems. Weber (1979) discussed the interchangeability of tandem/M/1 queues 
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in series. Boxma (1979) worked a tandem Queueing model with identical service times at 
both counties. Reiser (1979) did a queuing network analysis of computer co mmuni cation 
network with window flow control. Kuehn (1979) analysed networks approximately by de- 
composition. Strelen et.al (1979) analyzed Queueing networks with blocking using a new 
aggregation technique. Takahashi et.al (1980) considered the general system, which is a 
combination of tandem, split and merge configurations which is the most complicated to 
analyze, assuming that effective service times follows an exponential distribution, they 
developed a set of simultaneous nonlinear equations that must be solved to get performance 
measures. Ammar (1980) did not modeling and analysis of unreliable manufacturing 
assembly networks with finite storages Latouche et.al (1980) provided efficient algorithmic 
solutions to exponential tandwm queues with blocking. Bure (1980) discussed the 
dependence of sojourn times in tandem m/m/s queues. Seraphin (1981) made considerable 
efforts in developing and analyzing the interdependent communication networks composed 
of more than one transmission line connected in tandem. 

Ammar et.al (1981) gave equivalence relations in Queueing models of 
manufacturing. Shanti Kumar et.al (1981) provided open queueing network models of 
dynamic job shops internate. Ohno (1981) worked on a test problem that has a large number 
of states in order to evaluate algorithms for solving discounted semi-markov decision 
processes and compare several variants of there algorithms. Perros (1981) analyzed a 
symmetrical exponential open queue network with blocking and feed back. Sauer (1981) 
gave an approximate solution for Queueing networks with simultaneous resource 
possession. Pinedo (1982) worked on the optimal order of stations tandem queues. Wolff 



(1982) discussed tandem queues with dependent service in light traffic. Fakinas (1982) 
studied the generalized M/G/K blocking system with heterogeneous servers. 

Altiok (1982) gave an approximate analysis of exponential tandem queues with blocking 
Bouanaka (1982) discussed approximating queueing networks and facility planning. 
Gershwin et.al(1983) did modelling and analysis of three stage transfer line with unreliable 
machines and finite buffers. Whitt (1983) worked on the Queueing networks analyzer. 
Reiman (1984) discussed an open Queueing networks in heavy traffic. Coffman et.al (1984) 
used diffusion approximation for computer communication system. Jenq (1984) gave 
approximations for packetized voice in statistical multiplexer in process. Ganz O. (1984) 
provided an analysis of multiple bus synchronous and asynchronous communication 
system. Langaris et.al (1985) studied three stage tandem queue with blocking. Flores (1985) 
applied diffusion approximations for computer communication networks. Marshal et.al 
(1985) derived statistical measures for mixed date traffic on a local area network Yao 
(1985) did modelling for a class of state dependent routing inflexible manufacturing 
systems. Disney et.al (1985) reviewed queueing networks and gave survey of their random 
process Altiok et.al (1985) did approximate analysis of arbitrary configurations of open 
Queueing networks with blocking. Carmelila (1985) worked for improving mean delay in 
data communication networks by new combined strategies based in the soft principle. 
Pollok et.al (1985) gave an approximation method for tandem queues with blocking with 
exponential and general service distributions. Brandwajn et.al (1985) also gave an 
approximate method for tandem queues with blocking due to finite waiting room. Harrison 
(1985) worked on normalizing constants in Queueing networks. Bermanoded. et.al (1985) 
studied about optimal server location on a networks operating as an M/G/l queue. 
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Greenberg (1986) discussed queueing systems with returning customers and the 
order of tandem queues. Perros et.al ( 1986 ) gave an approximate analysis of open networks 
of queues with blocking having tandem configuration. Vinod et.al (1986) worked for 
approximating unreliable Queueing networks under the assumption of exponentially. 
Buzzacott et.al (1986) considered Queueing networks models of flexible manufacturing 
system. 

Altiok (1986) studied open networks of queues with blocking split and merge 
confiquration. Kerbache (1986) provided the generalized expansion method for open finite 
Queueing networks. Perros et.al (1986) discussed a computationally efficient 
approximation algorithms for analyzing open Queueing networks with blocking and having 
tandem configuration. Gershwin (1987) gave an efficient decomposition method for the 
approximate evaluation of tandem queues with finite storage space and blocking. 

Tsoucas et.al (1987) discussed the interchangability and stochastic ordering of /M/1 
queues in tandem. Shanti Kumar et.al (1987) derived bounding performance of Tandem 
queues with finite buffer spaces. Harrison et.al (1987) worked on Brownian models of 
open Queueing networks with homogeneous customer populations Altiok, et.al (1987) gave 
an approximate analysis of arbitrary configuration of open Queueing networks with 
blocking. Kerbache et.al (1987) provided the generalized expansion method for open finite 
Queueing networks. Schmidt (1987) found joint queue length characteristics in infinite 
server tandem queues with heavy traffic. 

Ohno (1987) computed optimal policies for controlled tandem Queueing systems. 
Brandwajn et.al (1988) gave an approximation method for tandem queues with bloking 
Vandijk (1988) obtained bounds for call congestion of finite single server exponential 
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tandem queues, Helber (1988) did decomposition of unreliable assembly, disassembly 
networks with limited buffer capacity and random processing times, Perros etal (1988) 
provided approximate analysis of product form type Queueing networks with blocking and 
dead lock. 

Bitran et.al (1988) discussed multiproduct Queueing networks with deterministic 
routing using decomposition approach and the notion of interference. Kerbache, etal 
(1988) studied asymptotic behaviour of the expansion method for open finite Queueing 
networks. 

King (1988) studied Queueing models for buffer with dial up servers. Gun, et.al 
(1989) provided an approximation method for general tandem Queueing systems subjects to 
blocking. Perros, etal (1989) gave a computationally efficient approximation algorithm for 
feed forward open Queueing networks with blocking. Akylildiz, etal (1989) found 
computational algorithms for networks of queues with rejection blocking. Harrison et.al 
(1989) did heavy traffic analysis of a simple open networks. Foss, (1989) derived some 
properties of open queuing networks. Brandwajn, et.al (1989) gave a note on approximate 
interactive solution of open tandem networks with blocking Lee, etal (1989) gave an 
approximate analysis for the merge configuration of an open Queueing networks with 
blocking Vandijk, et.al (1989) provided simple bound and monotonicity results for 
multiserver exponential tandem queues. Van Dijk (1989) obtained simple bounds and 
monotonicity results for multiserver exponential tandem queues. Jon and Perros (1990) 
provided an approximate analysis for open tandem Queueing networks with blocking and 
general service times. 
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Harrison et.al (1990) used the QNET method for two moment analysis of open 
queueing networks. Sigman (1990) discuss the stability of open queueing networks. Dallery 
(1990) gave an approximate analysis of generalized open Queueing networks with 
restricted capacity. Harry and Perros (1990) provided approximation algorithms for open 
Queueing networks with blocking. 

Christopher (1990) worked on diffusion approximations for computer 
communication networks. Gupta et.al (1990) considered production systems with two job 
classes and change over times and revisitation. Peterson (1991) worked on diffusion 
approximations for networks of queues with multiple customer types. 

Rybko, el.al (1991) discussed ergodicity of stochastic processes describing the 
operation of an open Queueing networks fouler and lelond (1991) studied local area 
network traffic characteristics with implications for broad band networks congestion 
management. Horvath (1992) gave strong approximations for open Queueing networks. 
Kella, et.al (1992) took a tandem fluid networks with levy input, in Queueing and related 
models. Bitran and Dasu (1992) reviewed open Queueing networks models of 
manuifacturing system. Dallery and Gershwin (1992) reviewed models and analytical 
results for manufacturing flow line systems. Jungshyrwiv (1992) did maximum entropy 
analysis of open Queueing networks with group arrival. Xinlichao, (1992) discussed 
reversibility of tandem queues with blocking. Iliadis (1992) studied synchronous versus 
asynchronous operation of a packet switch with combined input and output Queueing. 
Yakuo Hyshida (1993) did through put analysis of tandem type goback-NARQ, scheme for 
satellite communications. Dallery and Frain (1993) discussed decomposition methods for 
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tandem Queueing networks with blocking. Kelly and Laws (1993) considered dynamic 
routing in open Queueing networks. 

Harrison, et.al. (1993) gave a Brownian models of multi class Queueing networks. 
Spiros and Dnozis (1993) worked on entropy maximized Queueing networks with blocking 
and multiple job classes. Reyravi, Hassan (1993) derived limiteg distributions of the 
blocking propability for circuit switched networks. Park, et.al (1994) gave an approximate 
analysis of discrete time tandem networks with busty and correlated input traffic and 
customer loss. 

Waller (1994) provided a Queueing network models for field service support 
system. XINLTCHAO, (1994) considered a priority tandem queue with no intermediate 
buffer, Jain and Smith (1994) studied open finite Queueing networks with M/M/CK parallel 
servers. Yannopoul as and Alfa (1994) gave a simple and quick approximation algorithm 
for tandem split and merge Queueing networks with blocking. 

Sakasegawa and Kawashima (1994) studied about the equivalence of three types of 
blocking in non-markovian tandem queues. Balsamo, el.al (1994) did survey of product 
form Queueing networks with blocking and their equivalences. Frostand and Melamed 
(1994) did traffic modeling for telecommunications networks. Erramilli, et.al (1994) 
studied fractal traffic flows in high-speed communication networks. Avram et.al (1995) 
used optimal control approach for fluid models of sequencing problems in open Queueing 
networks. 

Dai (1995) discussed stability of open multi class Queueing networks via fluid 
modles. Nguyen (1995) used fluid and diffusion approximations of a two station mixed 
Queueing networks Pots, G. (1995) studied state dependent queuing networks. Liu and Lin 
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(1995) gave an efficient two-phase approximation method for exponential tandem 
Queueing systems with blocking. Gerasimov and Alxander (1995) worked on normalizing 
constants in multiclass Queueing networks. 

Vijaya Kumar (1996) studied some interdependent tandem Queueing models with 
applications to communications networks. Chakka and Mitrani (1996) gave approximation 
solution for open networks with break down and repairs. Gansnoah and Ryzinvangarrett, 
(1997) discussed about optimal control of a multi class flexible Queueing system. Lee, et.al 
(1998) worked on performance evaluation of open Queueing network with arbitrary 
configuration and finite buffer. Skianis, et.al (1998) studied arbitrary open Queueing 
networks with server vacation periods and blocking. 

Williams (1998) used diffusion approximations for open multiclass Queueing 
networks Dai, et.al (1999) proved heavy traffic limit theorem for a class of open Queueing 
networks with finite buffers. Knessi and Tier (1999) considered two tandem queues wit 
general renewal input using. Diffusion approximation and integral representation. Daduna 
and Szekli (2000) discussed the correlations of sojom times in open networks of 
exponential multiserver queues. Grassmann and Dekic (2000) gave an analytical solution 
for a tandem queue with blocking. Victar and Ivnitski (2001) considered networks of single 
server queues with dependent service times. 



1 .2 BROAD OUTLINES OF THE WORK 

Murty (1993) extended the single server interdependent queueing model to bulk 
service models M/M^Vl. He considered a different problem which is a generalization of the 
earlier Queueing model. Two variations are considered in the first one, the customer are 
served K at a time except, when less than K are in the system and ready for service at which 
time all units are served and in the second one the batch size must be exactly K and if it is 
not the server remains idle until the queue, size reaches K, a typical situation of this model 
is transportation of men and material are the customer. For this type of situation with bulk 
service rule they make use of along with the other assumption he assumed that batch size 
must the dependence structure given by Rao (1986). Further, Murty (1993) considered a 
slight variation of the model be exactly K and if not the server waits until such time to start. 
In order to reduce the idle time of the server it is feasible to consider the interdependence in 
the model also. He further assume that the number of arrivals and the service of any branch 
are correlated and follows a bivariate Poisson process. 

In the bulk service queueing models M/M M /1, Bailey (1954) and Jaiswal (1960) 
considered models in which the units arrive at random from a single queue in order of 
arrival and are used in batched, the size of each batch being either a fixed number of 
customers or the whole queue length which ever is smaller. Jaiswal (1961) extended this 
model to the case where at a service epoch m (0<M<S) if m persons are already present 
with the server then (S-m) persons or the whole queue length which is smaller will we taken 
in to service. The service rule is termed as Bailey’s is bulk service. However, in these 
models the arrivals and service process are poissonian and erlangian respectively with 
interdependent arrival and service processes. In both the models the system behaviour is 
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analyzed by obtaining the difference differential equation of the model and solving them 
through generating function techniques. The system characteristic like mean queue length, 
variability of the system size, coefficient of variations are derived and analyzed in the light 
of the dependence parameter. These models also include the earlier models as particular 
cases for specific values of the parameters. Murty (1993) considered the single server 
system with interdependent arrival and service process having the Bailey’s bulk service 
with variable capacity. Here he assumed that the server serves only at instants. 

ti,t 2 t n if m persons (0<M<S) are present in the waiting line at time t n 

then the server takes a batch of (S-m) persons or whole queue length which ever is smaller, 
where s is the service capacity. He also considered two server inter-dependent Queueing 
model with Baily’s bulk service rule for analysis. This model is an extension of the models 
considered earlier to server case it is assumed that the bulk service process having two 
service facilities with capacities bj and b 2 A batch of bi units or the whole queue length 
which ever is smaller is taken from the head of the queue for service in the first channel 
when even it is free. Similarly, the second channel on becoming free takes b 2 (<bi) or the 
whole queue length which ever is less. If both the server are idle and there is no queue the 
next unit to arrive always goes to the first service facility. Further he assumed that the two 
service facilities are independent of each other and the arrival and service process are 
independent. This sort of situation are common on marshalling yards with two engines, 
elevator process with two lifts etc. This model is also extended to multi server inter - 
dependent Queueing model with Bailey’s bulk service. In both the models he first develop 
the difference - differential equations and solve them through generating function 
technique. The system characteristics are derived analyzed in the presence of the dependent 
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parameter. The estimation of the parameters involved has not been consider which we 
propose to study in the present study and analyze the model with numerical idea. 

The work done on the discrete time transient solution for some Queueing models 
like multiserver queue with balking and reneging first passage time distribution, double 
ended queue, Busy period analysis queues with heterogeneous server and machine 
interference model etc. has been considered. 

Several Queueing problems have been solved using steady state conditions. As 
compared to those problems, it seems that not much have been done to obtain the 
corresponding transient solutions this is because of the fact that the transient solutions are 
not only mathematically intractable or excessively laborious but also computationally very 
costly. Therefore, we can say that most of Queueing theory results has concentrated on 
steady solution or some approximation. In most of the cases even steady state solution are 
difficult to compute. 

Chaudhary, Kapur et.al (1991,1992) have set a new trend in the numerical 
computations of models in queue through the technique of using roots. Closed form 
solution as well as exact computation results are obtained by this approach. Takac’s (1962) 
given two solution for M/M/l/oo neither of which is easy to compute. 

The first solution is in terms intervalas where as the second involves and infinite 
sum of Bessel frictions. The solution becomes a bit simpler if the waiting space is finite 
which may be true in many applications. In their case arrival and service rates are constant. 
Besides this they make use of spectral decomposition which require to find left and right 
eigen vectors. This is not easy if the matrix is large. The method developed by Chaudhry, 
Kapur et.al (1991) avoids spectral decomposition and well suited method for small and 
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large matrices. Numerical technique such as Runga- Kutta, Euler, Taylor and 
Randomisation have been used for finding transient solution- whereas the first there have 
been employed in solving differential equation the later on is particular suited for solving 
Queueing problems. In order to get greater accuracy one needs to increases to number of 
steps. Recently , Sharma and Das (1988) have obtained ansient solution to a special 
category of Markovian model using eigen values of matrices in Queueing theory. Standard 
numerical packages are used to obtain eigen values, which are difficult to obtain when the 
matrices are large and as a result computational difficulties arises. Chaudhary et.al (1991) 
have made attempt to obtain similar results is discrete time for finite waiting space 
problems in Queueing theory since the transient solution depend on the initial state of the 
system, it is interesting to know the effect on the system behaviour. Kobayshi (1983) has 
discussed the several systems which operate at discrete times of processor and several other 
examples in computer science. Chaudhry et.al (1991) give the transient solutions for 
general class of discrete times models in Queueing theory. In this way they assumed that 
the queue consists of finite waiting space. Inter arrival and service probabilities are 
dependent on the state of the system, the inter arrival and service time are geometric by 
independent of time and queue discipline is first in first out (FIFO). They obtained 
difference differential equations under these assumptions and solved these equations by the 
method of probability generating function. Cramer Rule has been applied for finding the 
solution of these equations. Explicit closed form expressions for distributions has been 
obtained in terms of the root of a characteristic equation. To find the eigen values or 
characteristic roots they make use of QROOT software package which is developed at 
Royal Military College at Canada by Chaudhry (1992). 
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Kapur, Sehgal el.al (1996) has obtained discrete time transient solution for the first 
passage time distribution under arbitrary initial conditions and finite waiting space. Further 
they worked on numerical computation of discrete time multi — server queue with balking 
and reneging. Although continuous time models are particular cases of discrete time 
models, this areas of research has remained neglected. Earlier attempts have largely 
concentrated on steady state solution of discrete time queues [Neuts (1989)]. It is in this 
sense that the study of discrete time models should stimulate to other areas such as 
computer science. Kapur, Sehgal et.al (1996) has further obtained discrete time transient 
solution for busy period analysis and double ended queuing system. They obtained closed 
form transient solution for the double ended Queueing system in discrete time. It is also 
further shown how the results in corresponding continuous time can be obtained further 
Kapur, Sehgal et. Al. (1996) has obtained a discrete time transient solution for Geom (N)/ 
Geom (N)/l/N machine interference model. 

It has been also discussed multi-operating system for the machine interference 
model with the and arbitrary initial number of machines interference. Extensive numerical 
results has also been given using the closed form expressions for multi - server machines 
interference model. The analytic and computational aspects of the model M [x VG/l and then- 
applications has been discussed by many author. The results have been discussed by many 
authors. The results have been reported in transform forms. The numerical aspects of 
Queueing models is at least as important as their analytical development. 

The work on discrete treatment solution for the first passage time multi server 
Queueing with bulking reneging busy period analysis has been obtained by Kapur P.K., 
Garg, R.B., Sehgal V.K. (1995) give numerical and transient solution of K out N:G systems 
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in discrete time with heterogeneous repair facilities. Kapur P.K., Sehgal V.K. , Garg R.B., 
Jha A.K. (1995) give numerical and transient solution of K ou,N:G system in discrete time 


with multi repair facilities. 



1.3-DISCRETE TIME TRANSIENT SOLUTION 

Several queuing problems have been solved using steady state conditions. As 
compared to these problems it seems that not much have been done to obtain the 
corresponding transient solution. This is because of the fact that transient solutions are not 
only mathematically intractable or excessively laborious but also computationally very 
costly. Therefore, we can say that most of queueing theory results has concentrated on 
steady state solution or some approximation. In most of the cases even steady state 
solution are difficult to compute. Chaudhry M. L., Kapur, P.K., and Templeton, 
J.G.C.(1991 ,1992) have set a new trend in the numerical computations of models in queue 
through the technique of using roots. Closed form solutions as well as exact computational 
results are obtained by this approach. 

Tackac’s (1 962)gives two solutions for the m/m/l/oo neither of which easy to 
compute. The first solutions is in terms of integrals whereas the second involes on infinite 
some of bessel functions. 

The solutions becomes a bit simpler if the waiting space is finite which may be 
true in any application. In their case arrival and service rates are constant besides this they 
make use of spectral decomposition which require to find left and right eigen vectors. This 
is not easy if the matrix is very large. 

The method developed by Chaudhry , M. L. Kapur , P. K. and Templeton , J. G. C. 
(1991) avoids spectral decomposition and well suited method for small and large matrices. 
Numerical techniques such as Ranga Kutta, Euler, Taylor and Randomization have been 
used to find transient solutions. Whereas first three have been employed in solving 
differential equation, the later one is particular suited for solving queueing problems. In 
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order to get greater accuracy one need to increase the number of steps. This together the 
number of simultaneous equations to be solved shows down the solutions process 
considerably Recently Sharma and Das (1988) have obtained transient solutions to a special 
category of Markovian model using eigen values of matrices in queueing theory. Standard 
numerical packages are used to obtained values, which are difficult to obtained when the 
matrices are large and a results computational difficulties arises. 

Chaudhry, M.L., Kapur , P.K. and Templeton, J.G.C. (1991) have made attempt to 
obtained similar results in discrete time far finite waiting space problems in queueing 
theory. Since the transient solution depend on the initial state of system, it is intersecting to 
know the effect on the system behaviour. Kabayashi (1988) has discussed the several 
system which operate at discrete time many machine cycling of a processar and several 
other examples in computer science. 

Chaudhry, M.L., Kapur. P.K. and Templeton, J.G.C. (1991) give the transient 
solutions for a general class discrete time models in queueing theory. In this they have 
assumed that the queue consists of finite waiting space. Interarrival and service 
probabilities are dependent on the state of the system, the inter arrival and service 
probabilities are dependent on the state of the system, the inter arrival and service time 
distributions are geometric but independent of time and queue discipline is first in first out. 

By using these assumptions firstly they make the difference equations. To salve 
these equations they use the method of probability generating functions write these 
equations in matrix form. Cramer rule has been applied for finding the solutions of these 
equations. Explicit closed form expressions for distributions has been obtained in terms of 
the roots of a characteristic equation. 
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To find the eigen values or characteristics roots they make use of QROOT software 
package which is developed at Royal Military College at Canada by M.L. Chaudhry (1992). 
For the analysis of the model the following notations are used 
Xk - Number of customers in queue at epoch k. 

N - Size of waiting space. 

X n - Inter arrival probability when n customers are in the system 
Pn - Service probability when n customers are in the system. 

— ftn (1 “ ^-n) 

Vn ~ (1 " ft-n) 

Here P m (n) denote the probability that the system is in the n th state at the beginning of m* 

epoch Xk, k >0 in an interger valued discrete stochastic process taking value (0,1,2 ,N). 

Xk=n (0 < n < N) implies that there are n customers in the system at epoch k. The difference 


equations are 

Pm+l (0) - P m (0) = - VJ/o P m (0) 4>! . ...(U) 

Pm +1 (n) - Pm (n) = - Pm(n) (<t>n + ^n) + Pm (n-1) vj/ n -i + P m (n+1) <j) n +l 

; 1 < n < (N-1) (1.2) 

Pm+l (N) — P m (N) ~ ~ cf>N Pm (N) + Pm (N-1) V}/ N -i (1.3) 


and Po (i) = i ; 0 < i <N 

Let P z (n) be the p.g.f. of P m (n) defined as 

OO 

P z (n) = E Z m Pm (n) !z|< 1 

m=0 

Now taking the p.g.f. of equation (1.1), (1.2) and (1.3), we have 
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Pm*l (0) z» - P m (0) z“ - - y, P. (0) z"+ ^ 

” ” ® » 

Z Pm+i (0) Z m -Z Pm(0) Z m = - Xj/oZ Pm (0) Z m + fa Z}p m (l) Z m 

m=0 m=0 m=n ^ ; ' 



(P z (0) - Po(0)) - p z (0) = - Vj/o P z (0) + (J), p 2 (1) 
- Pz (°) ' -L p o(0) = - \j/ 0 P z (0) + cj), p 2 (i) 


(s+Vo) Pz (0) - 4> P z (1) = _L_ P 0 (0) 

Wher s = 1 -z 

z 

00 CO oc 

Z Pm+i (n) z -Z Pm(n) Z m = Z Pm (n) z m (-(j) n -V|/ n ) + vpn.! Z P m (n-l)z 
m_0 m= 0 m=0 m=0 

OO 

+ <f>n + l Z Pm (n-1) Z m 
m=0 

~z ~ {Pz(n) ~ Po(n)} ~ Pz ( n ) = C-<f»n -V|/ n ) P z (n) + Vj/n.! P z (n-1) +<J) n+ , p z (n+l) 


. Y~] Pz (n) ’ ° (n) = ^ n ‘ ¥n) Pz (n) + ¥n ' 1 Pz ( n_1 ) + <f>n+i Pz (n+l) 

(s +(()„ +y n ) P z (n) - Vn ., P z (n-1) - (j) n+1 p z (n+l) = (l/ z ) P 0 ( n ) 

-Vn-I Pz (n-1) + (s+<j>„ + Ij/ n ) P z (n) - (j) n+1 p 2 (n+l) = (i/ z ) p 0 ( n ) 

; 1 <n^N- 1 

00 00 °0 oo 

Z P m+ 1 (N) Z m -Z Pm(N) Z m =- (j) N Z Pm (N) Z m + ZPm (N-l) Z m . 
m_0 nx=0 m=0 m=0 

( l/z){P z (N) - P 0 (N)} - P z (N) = -<|> N P Z (N) +*P n .iP z (N-1) 


r ^ 
1-z 
z 


Pz(N) (l/z)P 0 (N) - - (j) N P Z (N) + \|/n_ i P z (N-1) 


( s+ ^n) P z (N) -\]/ N -iP 2 (N-1) = 1/z P 0 (N) 

- v)/ N-i Pz (N-l) + (s + 4> n ) P z (N) = (1/z) P 0 (N) 


Where s = i-z 
z 


( 1 - 6 ) 


Now we have 
A.P. = [Sko,§ici,- 


W 


.(1.7) 


Where A is real tridiagonal (N+l) * (N+l) Matrix, p is colunm veetor and S ki is a 
kronecker delta defined as 


Ski = 

Assuming po = 0 and X N = 0 

0 1 2 


1/z 

0 


k = i 

otherwise 


A(s) : 


N-2 N-l 


N 


0 

s+yo 

H>i 

0 

- - - 0 o 

0 

1 

-Vo 

s+^i+Vi 

- <j>2 

- - 0 o 

0 

N-l 

0 

0 

0 

■ -VN-2 S+ljjN.j + 

-<|>N 

N 

0 

0 

0 

■ ■ ■ -Vn-1 

S+(|) N 

and 





/I OY 


Pz(0) 

Pz(l) 


L p z (N)J 


(N+l) x 1 


From (1.8) Cramer’s Rule P z (n) are explicitly determined as 


I A n +i (s) | 

P z(fl) = ' ; I z | < 1 

I A(s) | 


(1.9) 


( 1 . 10 ) 


Applying some row and column transformation on I A(s) | it may be expressed as 
s I D(s) | , where D(s) is the matrix of order NxN given by 


1 2 3 - - > N-2 N-l N 

1 l+<j)i+H/ 0 W(<j>iH/i) 0 .... - 

I W(<J)i l|/i) S+(j) 2 +V|/i -V((|)2 vj/ 2 ) - - - - - - 


N-l 


- - W(<j>N-2 VN-2) 


S+(^N-2+ VN-2 


V(())n-1 VN-i) 


N 



W(<j)N-l VN-l) S+<j)N+ VN-1 


If D(s) is expanded it will be a polynomial of degree N. The roots of the polynomial 
I A(s) | are real, negative and distinct (one root being zero). Let at (k = 0,1,2.....N) be the 
root of | A(s) | with 


ao = 0 then 


and hence 


N 

I A(s) | = s n (s - at) 
k= 1 


I A n+ i(s) | 

p z( n ) = ~ 0< n <N (1.11) 

S FI (S-CCk) 

k= 1 

resolving the right hand side of P z (n) into partial fraction replacing s by (l-z)/z using initial 
condition and comparing coefficient of z m , we have 


i-1 N 

^ FI <|> rfi I a kr ( 1+ a k ) m 
r=n k=l 


P m (n) = b n + I 


N 

E a kn ( 1+ a k f 
k=l 


n-1 N 

FI vp r Ea kn (l+a k ) m 

r=i k=l 


where a kn ’ S and b n ’ S are defined as 

Cn-i (oc k ) D n (cc k ) 
N 

a k n (a k - <Xj). 
j=l 

C N -i (a k ) Dj (a k ) 


tX-kn 



N 

a k n (otk-aj) 

j=l 


C N -n (ak) Pi (ak) 
N 

a k n (ak-ccj) 

j=l 


0 < n < i 



i<n<N 

. ... ( 1 . 12 ) 


0 < n < i 


n = i 


i <n<N 


and 


u = CN-n (0) D n (0) 

D n — 

N 0 <n<N 

FI ( - a k ) 

k=l 

(l- 14 ) 

where C n (s) and D n (s) being the determinants obtained by the bottom right and top left and 
(n x n ) square matrices formed from A(s) such that 
I A(s) | = Cn+i (s) = Djm+i (s) 

for convenince, we write C n (s) and D n (s) as C n and D n respectively then, 

Ci= [s+<j)N+l-i + V(/N+l-i] Cj.i - [(j) N+ 2-i VN+l-i] Cj-2, 2 < i <N+1 (1.15) 

And 


Dj = [ s +<j) i.j+ v[/ ,.]] Di_i - [<)> i-rH/ i- 2 ] Dj .2 2<i <N+1 (1.16) 

Using QROOT software packages, we find the root a k called the characteristic equation of 
A(s) after finding the roots they discuss many cases and find the numerical results. 
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1.4- DISCRETE TIME TRANSIENT SOLUTION FOR ■ 
GEOM (n) /GEOM (n) 12 /N WITH HETEROGENEOUS SERVER 

We assume that the inter arrival probabilities and service time probabilities of first 
and second server to be geometrically distributed with parameter X, pi, and p 2 respectively. 
We also assume that pi<p 2 that is the service time probability for first server is less than 
that of second server, we are considering modified queue discipline i.e. the first arriving 
unit from amongst the initial number of unit present of the start of the service join the first 
counter for service. 

Therefore the arriving unit goes to the counter which it find free, the maximum 
number of customer in the system is restricted to N. We further assume that there is no unit 
initial waiting at the time t=o when the service starts 
Notation:- 

Xk> denote the number of customer at epoch k 
N:- Size of waiting space. 

X:- interarrival probability of a customer 
pp- service probability of a customer for server one 
p 2 :- service probability of a customer for server two 
<j>i= pi(l-X) ,<|>2 = P2(1-^), 'F=X(l-pi-p 2 ) 

Analysis of the model : 

Let P m (n) (n=0, 1, 2; ~N ) denote the probability that the system is in the n* state at 
the beginning of the m th epoch or time slot. Let X k be the number of customers in the 
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system at discrete time epoch k. Then X k, k >0 is an integer valued discrete stochastic 
process taking values 0,1, N. Xk = n (0 < n< N) implies that there are n 


customers in the system at epoch k . 

The following difference equations can be written as 

Pm+l (0) - P m (0) ^ Pm(0) + <|)i P m (1) ....... (1.17) 

Pm+l (1) -Pm (1) = - Pm (1) (y + 4»1 + h) + k P m (0) + ( «j>l +(j) 2 ) P m (2) (1.18) 

Pm+l (n) -Pm (n) = - Pm (n) (y + 4>1 + fe) + y P m (n- 1 ) + ( <j>i +<j) 2 ) P m (n+ 1 ) (1.19) 

2 < n <N-2 


Pm+l (N-l) -Pm (N-l) = - P m (N-l) (q/ + 4>1 + (j) 2 ) + q/ Pm (N-2) + (p, + p 2 ) P m (N) -(1-20) 

Pm+l (N)-P m (N) = q/Pm(N-l)- (P! + p 2 ) Pm(N) (1.21) 

Where Po(i) =1 ; 0 < i < N 

Let P(n) be the steady state distribution i.e 

Lim P m (n) = P(n) 
m^o) 

Let P z (n) be the probability generating function (p.g.f.) of P m (n) defined as 

CO 

G(z,n) = P z (n) = X z m Pm (n) I z | < 1 

m=0 

Taking the p.g.f of equation (1.17) to (1.21) form this we multiply the equation by z m and 
taking summation from 0 to oo for m and using (1-z) /z = s, we get 


(s+A.) P z (0) - <j)i P z (1) = 1/z (1.22) 

-X P z (o) + (s+m/+ h+h) Pz (1) - (<t>i^ 2 ) P 2 (2) = 1/z - (1 .23) . 

-q/ P z (n-l) + (s + H/ +<)>i+<j> 2 ) P z (n) - (4>i+<t>2) Pz(n+1) = 1/z (1-24) 

2< n <N-2 
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-y P z (N-l) + (s + f + 4> 2 ) Pz(N-l) - (m+ p 2 ) P z (N) = 1/z .........(1.25) 

P 2 (N-l) + (s + P 1 +P 2 ) P z (N) = 1/z (1.26) 

These equation (1 .22) to (1.26) can be written as the matrix form 

AP = [8 k0 ,5 k i, SkN]' (1.27) 


Where A is a real tridiagonal (N+l) x (N+l) matrix , P is a column vector of order 
(N+ 1) xl and 8 k j is the Kronecker Delta defined as 
P 1 /z ; 



0 ; otherwise 


4 


We have 

0 

1 


0 1 2 

s+X -(j>i 0 

-X S+\|/+(j)i+<j)2 - (4>1+<|)2) 

0 -y (s+<t>i+(j) 2 ) 


A(s) = 


- N-2 N-l N 

- 0 0 0 

- 0 0 0 


N-2 - 

N-l 0 0 0 

N O 0 0 


- - - (s+Vl/+(|)i+(j)2) -(4>1 +< I>2) ■ 0 

0 S+V|/+(|)i+<j)2 -(|Ii+q.2) 

- - - 0 -\j; S+m+|t2 


And P = [P z (0), P z (1), Pz (N)] ; 

From equation (1.27) , using Cramer’s rule P z (n)are explicitly determined as 

Pz (n) = 1 A n +i (s) 1 ; 0< n <N 

I A (s) | 
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D(s)= 


Where A n +i (s) is obtained from A(s) by replacing the (n+l) th column of A(s) by the 
right hand side of equation (1.27) and I A (s) | is the determinate ofA(s) 

Applying some row and column transformation on | A (s) | , it may be expressed : 
s | D(s) | is a real symmetric tridiagonal matrix of order (NxN). 
s+(<j>i+X) -V(<j>i \| /) 0 - 0 0 0 

-V(c|)i\| /) s+((j)i +([)2 )-H|/ -Vv|/ - - - 0 0 0 

Oh+fc) 


Wvj/ (<j)i+(j) 2 ) s + l P+((j)i+(j) 2 ) 




S+'F+(|J,i+(4.2) 


I D (s) | is a polynomial of degree N of s. The roots of | D(s) | are the negatives of the 
eigenvalues of the matrix D(0). Hence the roots of the polynomial | A (s) | are real, 
negative and distinct . Let (Xk (k — 0, 1,2, N) be the roots of | A(s) | with ao=0 then, 


I A (s) | = s n (S - a k ) 
k= 1 


and hence 


P z (n) 


I A n +l (s) | 
N 

s FI (s-aj) 

J = 1 


0<n<N 
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Resolving the R.H.S of P z (n) in to partial fractions, replacing s by (l-z)/z using 
initial conditions and comparing coefficient of z m , we can find out the value P m (n). It will 
be bounded for 1 1+a k |<1. 

This work in that sense gives impetus to the analysis of discrete time models with 
heterogeneous servers. 



1.5 DISCRETE-TIME TRANSIENT SOLUTION FOR A FIRST 
PASSAGE TIME DISTRIBUTION IN QUEUEING THEORY 

Most of the queueing theory literature concentrates on finding the steady state 
solution or approximations. Chaudhry M.L., Agarwal, M. and Templetion, J.G.C. (1992) 
have mostly concentrated on this using the technique of roots. Earlier attempts at finding 
transient solution can be attributed to Takacs, L (1962) and Morse P.M. (1958). However, 
there are computational difficulties with their methods. 

Now with the increased skill available in computations with the use of computers, 
researchers especially in computer science have started looking for transient solutions and 
easy to compute closed form solutions. Recently, Sharma, O.P. and Das, S. (1988) have 
provided transient solutions to a class of Markovian models in Queueing Theory. 

The transient solution is not independent of the initial state of the system it is 
interesting to know its effect on the system’s behaviour. Some system may not exist long 
enough to reach their steady state. 

There are several systems which operate at discrete times see Kobayashi H. (1983) 
As a result, it becomes important to study them. 

A discrete time Queueing model for the first passage time distribution to a 
absorbing state given the initial state. We give closed form solution to this class of 
problems in terms of the root of a polynomial in z-transform and results are computed even 
when the matrix involved are large. It is also shown, how the results for the continuous case 
can be obtained. 

Result presented in this section further unify the treatment given by Chaudhry M.L., 
Kapur P.K., Templeton J.G.C. (1991). 



We assume that 


1. The queue consists of finite waiting space. 

2. Inter-arrival and service probabilities are dependent on the state of the system. 

3. Inter-arrival and service time distributions are geometric but independent of 
time. 

4. Queue discipline is first in first out. 

The notations used are : 

Xk : Number of customer in the system at epoch k. 

N : Size of the queue (maximum) 

X n : Inter-arrival probability when n customer are in the system. 

p n : Service probabilities when n customer are in the system. 

Vn : n (1 - Jin) 

A 

<j) n - n (1 _ An) 

h : Absorbing barrier (h >0) (h<n) 

ANALYSIS OF THE MODEL : 

Let Xk be the number of customers in the system at discrete-time epoch k. 
Then Xk, k>0 is an integer-valued discrete stochastic process taking values 

(h,h+l, N). Xk = n (h<n<N) implies that there are n customers in the system at 

discrete time epoch k. When a customer arrives or leaves, a discontinuity in the stochastic 
process occurs. The process Xk behaves as a discrete-time Markov process and represents 
^ the state of the system. 
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Denote the probability that the system in the state n at the beginning of the m th 
epoch as P m (n), (n=h,h+l ,N) 

Backward First Time-Distribution Analysis : 

Consider the case for the elapsed time slot m, the following difference equations 


may be easily written before the queue length for the first time reaches h. 

Pnv+1 (h) — Pm (h) = <t>h+l Pm(h+1) .....(1.28) 

Pm+1 (h+1) - P m (h+1) = - (Vl/ n +l+<t>n+l) Pm(h+1) + <l>h+2 Pm (h+2) .(1 .29) 

Pm+l (ft) — ?m (h) ' ('f n"Hv) P m(h) + l|ln-l P m (H-l) + <j>n+l Pm (h+1) .(1 .30) 

(h+2)<n< (N-l) 

Pm+l (N)- Pm (N) = - hi Pm(N) + Vn ., P m (N-l) (1.31) 

Where A,n=0 and Po(i) = 1 , h < i < N 

Let P n be the steady state distribution, i.e. 

Lim P m (n) = P(n). 

m->oo 


If such a distribution exists, it is unique. Solving (1.28) to (1.31) for the stationary 
case, we get . 

P(h) = 1, P(1) = 0, h+1 < 1 < N 

Let Pz(n) be the probability generating function (p.g.f.) of P m (n) defined as 
oo 

G(z,n) = P z (n) = Z z m P m (n), I z I < 1 

m=0 

Taking the p.g.f. of equations (1 .28) to (1 .31), we have 

AP = [8kh, 8ic(h+i) 8kN] / (1.32) 

Where A is real tri-diagonal (N-h+1) x (N-h+1) matrix, P is a column vector and Ski 
is the kronecker delta defined as 



Ski - "s 


1/z 

0 


k=i 

otherwise 


Definning s= (l-z)/z, we have 
h h+1 h+1 


h+3 


h 

h+1 

h+2 


A(s) : 


N-l 

N 


s -<|>h+i 0 

0 S+VJ/ h +l+<t>h+l -<(>h+2 

"Vh+1 S+V(/hH-2+(|)h+2 -<t>h+3 


P = 


Pz(h) 

P z (h+1) 


Pz(N) 


(N-h+1) x 1 


N-l 


N 


VJ/N-2 S+\|/n-1+<|>N-1 


"W-l S+([>n 

(N-h+1) x (N-h+ 


) 


From equation (1.32), using Cramer’s Rule P z (n) are explicity determined as 


A n .h+i(s) 


Pz(n) 


I A(s) 


h<n<N 
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Where A n -h+i is obtained from A(s) by replacing the (n-h+l) th column of A(s) by 
R.H.S of (1.32) and | A(s) | is the determinant of A(s). 


Applying some row and column transformation on | A(s) 

I , it may be expressed as 

S 1 D(s) | , where D(s) is a real, symmetric, tri-diagonal matrix 

of order (N-h) 

x (N-h) 

specifically D(s) = 



s+(<t>h+0 W(^ IVW ) - - - - - 


0 

W(<j)h+lM'h+l) S+((f) h + 1 +V|/h+ 1 ) - - - - - 

- ■ 

0 

V(<ktWh+2) 

- 

- 

■ ' • -V(<t>N-2 W- 2 ) 

S+C^N-l+W-l) 

-^N-l^N-l) 

- - - - 

i 

© 

- I 

s+^n+Vn-i) 


D(s) is a polynomial of degree (N-h) in s. It may be noted that the roots of | D(s) | are the 
negatives of the eigen values of the matrix D(0). 


It may be observed that D(0) is a positive definite, symmetric, tri-diagonal matrix. 
The roots of the polynomial |A(s) | are real, negative and distinct (one root is zero). 

Let ak (k=0,l,2 ,N-h) be the roots of | A(s) | with ao = 0. Then. 


h<n<N 
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N-h 

|A(s) | = s II (s - ak) 
k=T 

And hence 

I A n -h+l(s) I 

Pz(n) = — 

N-h 

s n (s - aj) 

j=i 



Resolving the R.H.S of P 2 (n) into partial fractions and replacing s by (1 - z)/z, using 
initial conditions and comparing the coefficients of z m , we get. 
i-1 N-h 

Pm(n) = 1+11*,, Z Skh (l+OCk) m 

r=hk=l 

i-1 N-h 

Pm(n) = n <t>„ Z a kn (l+cc k ) m , h < n < i 

r = n k = 1 

N-h 

Pm(n) = Z a kn (l+a k ) m , n=l 

k= 1 

n-1 N-h 

Pm(n) = n Wr+1 Z a kn (l+a k ) m , i<n<N 

r = i k= 1 


h < n < 1 


n=l 


C N .„ (a k ) Dj-h (a k ) 

a kn — 

N-h 

a k n (oc k - aj) 

j=l^k 


' ■■ .‘1 

" ' 1 

„• 1 ' i 


CN-i (cX-k) 



N-h 



a k 

n 

(a k - 

aj) 

j : 

=l*k 



Cn- 

i (ak) 

D n -h 

(ak) 


N-h 



a k 

n 

(a k - 

aj) 

j : 

r-H 

II 



Cn- 

i(a k ) 

Dn-h 

(ak) 


N-h 



a k 

n 

(a k - 

aj) 

j 

=l*k 
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With C n (s) and D„(s) being the determinants obtained by the bottom right and top 
left (nxn) square matrices formed from A(s) such that 
I A(s) | = C]\i-h+i(s) = D N . h+1 (s) 

C n (s) and D n (s) may be determined by the following recurrence relations. 

Assuming Co(s) and Do(s) = 1 , Ci(s) = s + ^ Di(s) = s and An = y n = ((>„ = 0 
Ci(s) = (s + \|/N+l-i + ^N+l-i) Cm - \j/N+2-i ^N+2-i Q-2 

; 2 < i <N-h+l 

Di(s) = (s + vph+i-l + <t>h+i-l) Di-l - v h+i-l <|> h+i -1 Di -2 

; 2 < i <N-h+l 

Using the standard IMSL packages one can find the eigen values and hence the 
zeros of the polynomial | A(s) | . Further using the recurrence relations for Q(s) and Di(s) 
can easily be evaluated. 

Since (l+ak) m ->0 as m->oo the steady state distribution P(n) is given by 
P(h) = 1. P(n) = 0 h+1 < n < N 

Important Performance Measures : 

Using closed form expressions for P m (n), some important measures can be 
analytically and numberically derived 

1 . Expected number of customers in the system (for fixed i) 

N 

E (X m ) = E n P m (n) 
n=h 

2. If Y m denotes the number of customers present in the queue. 

N 

E(X m )= S (n-r)P m (n) 
n=h+r 
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3. Probability that the system state is greater than a given number c is given by 


N 


(c>h) I P m (n) 


n=c 


4. Relaxation time which is a measure of length of time. 


RT = 1/1 
If m»RT, P m (n) = P(n) 


r N ^ 
min 
i=l 


(-Re(aO) 

Vn 


Forward First Passage Time : 

Next we consider the absorbing barrier on the maximum queue size, we may write 
the difference equations as 
Pm+l (0) - P m (0) = (v}/ 0 ) P m (0) 

Pm+i(n) - P m (n) =■ - (V|/ n +<t>n) Pm(n) +\\i n -i Pm(n-l) + <t>n+i P m (n+1) 1< n < N-2 

Pm+l(N-l) — P m (N-l) (V)/]\j.i+<j)!v(.i) P m (N-l) +Vj/)vj.2 Pm(N*2) P m +i(N) 

Pm+l(N) — Pm (N) tyN-1 Pm(N-l) 

Where }io = 0 and Po(i) = 1, ;0<i<N 

Proceeding as above the steady state probabilities may be given as 
P(i) = 0 ; 0 < i <N-1 P(N) = 1 

The probabilities P m (n) may be expressed as 
i-1 N 

P m(n) = 1 1 4*r+l &kn (l"N^k) ’ 0 < n < i 

r=n k= 1 

N 

Pm(fr) = 3kn(l"^0tk) » j tt — i 

k= 1 
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n-1 N 

Pm(n) = n Vr S akn(l+ock) m , ; i <n<N 

r = i k = 1 

N-1 N 

Pm(N) = 1 + n Vr 2kN (l+CCk)” 1 
r = i k= 1 

Cn-j (<Xk) D n (otk) 

akn = . — , ; 0 < n < i 

N 

ak n (ak-oj) 

j=ia 


C N -i (®k) Dn (&k) 

a kn = . ; n =i 

N 

ot k n (ok - oj) 

i=i*k 

C>j-n (otk) Dj (otk) 

a k n = , ;i<n<N 

N 

ak n (a k - «j) 

j=l*k 

Di(a k ) 

a kN = • , 

N 

a k n (otk-aj) 

j=l*k 

C n (s) and D n (s) are defined earlier, 
i- CHANNEL BUSY PERIOD : 

We define i-channel busy period (o < i < N) to begin with an arrival to the system at 
an epoch when there are (i-1) customers in the system to the very next epoch when there are 
again (i - 1) customers in the system. Assuming X n and p, n to be the inter-arrival and service 
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probabilities respectively. When there are n customers in the system, the following 
difference equations may be written 

Pm + 1 (i-1) - P m (i-1) = 

Pm+l(i) - Pm(i) = - (Wi+<t>i) Pm® +<t>i + 1 P m (i+1) 

Pm+i(n) - Pm(n) = - (M/n+<t>n) P m (n) +\|/ n -i P m (n-1) + (j) n+ i P m (n+1) i < n < N-l 
Pm + l(N) - P m (N) = - <j) N P m (N) + l|/N-l Pm(N-l) 

Where Xn = 0 and Po(i) = 1 

The solution of these equations can be obtained as before with h=i-l and.h+l=i. 
Pm(i-l) and <j>i P m (i) are the probability distribution and probability mass function of the 
busy period. 

NUMERICAL RESULTS : 

The numerical results for only moderate values of N are given though were no 
problems even for large values of N. 

Case (i) Backward First Passage Time : 

Discrete Case : Assume r = 6, N=20, X=0.8, p=0.15, m=T0, i = 4,5,6, ...,20, h=4 

(h < i < N) and Po(i) = 1/17. Table 1.51 gives the probabilities P m (n) for different i, the 
unconditional probabilities Q m (n) and steady state probability P(n). The last two rows give 
the values of E(X m ) and E(Y m ). For i=20, m=700 which is »RT=47. 

Continuous Case : Assume r=6, N=20 X,=0.8, p=0.15, t=10, i=4,5,6,......,20. h=4 

(h < i < N) and Po(i) = 1/17. Table 1.52 gives the probabilities P n (t) for different i, Q n (t) and 
P(n). For i=20 the time t=900 which is »RT = 86. 
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Case (ii) Forward First Passage Time : 

Discrete Case : Assume r = 6, N=16, X=0.8, p=0.15, m=10, i=0, 1,2,... ...... ,16, and P 0 (i) - 

1/17. Table 1.53 gives the probabilities P m (n) for different i,Q m (n) and P(n). 

Continuous Case : Assume r=6, N=16 X=0.8, p=0.15, t=10, i=0,l,2. ,16. 

Po(i) - 1/17. Table 1.54 gives the probabilities P n (t) for different i, Q n (t) and P(n). 

Case (iii) i-channel busy period : Assume r = 6, N=20, A,=0.8, p=0.15, m=10, i=5,6. 
Table 1.55 gives the probabilities P m (n) and the steady state probabilities. 




Table 1.51 : Probabilities P m (n) P(n) and means for Geom(n)/ Geom(n)/r/N with r = 6 

N= 20, k= 0.8 (r=0.15 m=10 i=4,5,6 ........20, h =4,(h <i <n) and Q m (n) with 

PoCO =1/17 


i/n 

4 

5 

6 

18 

19 

20 

Qm(n) 

P(n) 

4 

1.0000 

1211111 '; 

0.3163 

0.0000 


0.0000 

0.1217 

1.0000 

5 

0.0000 

0.1367 

0.1874 

0.0000 


0.0000 

0.0370 

0.0000 

6 

0.0000 



0.0000 


0.0000 

0.0717 

0.0000 

7 

0.0000 

0.0789 

0.1376 

0.0000 


0.0000 

0.0627 

0.0000 

8 

0.0000 


0.0444 

0.0000 


0.0000 

0.0597 

0.0000 . 

9 

0.0000 

0.0036 


0.0000 


0.0000 

0.0590 

0.0000 

10 

0.0000 

0.0005 

0.0017 

0.0000 


0.0000 

0.0589 

0.0000 

11 

0.0000 

0.0000 

0.0002 

0.0003 

0.0000 

0.0000 


0.0000 

12 

0.0000 

0.0000 

0.0000 

0.0022 

0.0003 

0.0000 


0.0000 

13 

0.0000 

0.0000 

0.0000 

0.0115 

0.0022 

0.0003 


0.0000 

14 

0.0000 

0.0000 

0.0000 

0.0421 

0.0115 

0.0023 


0.0000 

15 

0.0000 

0.0000 

0.0000 


0.0421 

0.0124 


0.0000 

16 

0.0000 

0.0000 

0.0000 


0.1103 

0.0466 

0.0583 

0.0000 

17 

0.0000 

0.0000 

0.0000 

0.2551 

0.2039 

0.1254 


0.0000 

18 

0.0000 

0.0000 

0.0000 

0.2134 

0.2619 



0.0000 

19 

0.0000 

0.0000 

0.0000 

0.1164 


0.3124 



20 

0.0000 

0.0000 

0.0000 




0.0217 

0.0000 

Sum 



1.0000 





1.0000 

E ( xm ) 


4.8954 

5.4449 


17.8967 

18.5685 

11.1196 

4.0000 

E ( Ym ) 





11.9867 

12.5685 

5.3999 

0.0000 

















Table 1.52 : probabilities P n (t) .P(n) and 


|i =0. 1 5, t =10, i =4,5,6, ..,20. h =4 (h 


i/n 

4 

5 

6 

18 

4 

1.0000 

0.8219 

0.6693 


5 

0.0000 

0.0153 

0.0269 


6 

0.0000 

0.0239 

0.0422 


7 

0.0000 

0.0279 

0.0489 


8 

0.0000 

0.0277 



9 

0.0000 

0.0244 

0.0452 



0.0000 

0.0179 

0.0371 


ii 

0.0000 

0.0145 

0.0282 


12 

0.0000 


0.0199 


13 

0.0000 


0.0131 


14 

0.0000 


0.0081 


15 

0.0000 


0.0048 

0.1995 

16 

0.0000 


0.0026 

0.1137 

17 

0.0000 


0.0014 

0.1236 

18 

0.0000 


0.0007 

0.1279 

19 

0.0000 





0.0000 



0.1181 

Sum 

1.0000 




E(X t ) 

4.0000 

4.8287 

5.5862 

15.964! 

E(Y t ) 


0.4879 

0.9517 

9.9698 


5 for M/M/R7N with r = 6, N =20,X=0.8, 
<i <n) and Q n (t) with p 0 (i) = 1/17, 

55 \m 


19 


20 















Table 1.53 probabilities p m (n) and means for Geom/Geom(n)/r/N with r=6 7,=0.8, 
p=0.15, m =10, i —1,2, 16, and Q m (n) with p 0 (i) =1/17. 


i/n 

0 

1 

2 

14 

15 

16 

Qm(n) 

p (n) 

0 

0.0004 

0.0003 

0.0002 

0.0000 

0.0000 

0.0000 

0.0001 

0.0000 

1 

0.0075 

0.0015 

0.0033 

0.0000 

0.0000 

0.0000 

0.0012 

0.0000 

2 

0.0510 

0.0377 

0.0269 

0.0000 

0.0000 

0.0000 

0.0093 

0.0000 

3 

0.1763 

0.1440 

0.1144 

0.0000 

0.0000 

0.0000 

0.0399 

0.0000 

4 

0.3255 

0.3006 

0.2707 

0.0000 

0.0000 

0.0000 

0.1014 

0.0000 

5 

0.3047 

0.3290 

0.0.3448 

0.0000 

0.0000 

0.0000 

0.1567 

0.0000 

6 

0.1219 

0.1611 

0.2034 

0.0000 

0.0000 

0.0000 

0.0427 

0.0000 

7 

0.0120 

0.0204 

0.0324 

0.0003 

0.0000 

0.0000 

0.0753 

0.0000 

8 

0.0007 

0.0017 

0.0036 

0.0022 

0.0003 

0.0000 

0.0614 

0.0000 

9 

0.0000 

0.0001 

0.0003 

0.0115 

0.0022 

0.0000 

0.0591 

0.0000 

10 

0.0000 

0.0000 

0.0000 

0.0421 

0.0144 

0.0000 

! 

0.0586 

0.0000 

11 

0.0000 

0.0000 

0.0000 

0.1097 

0.0411 

0.0000 

0.0579 

0.0000 

12 

0.0000 

0.0000 

0.0000 

0.2015 

0.1046 

0.0000 

0.0552 

0.0000 

13 

0.0000 

0.0000 

0.0000 

0.2527 

0.1832 

0.0000 

0.0478 

0.0000 

14 

0.0000 

0.0000 

0.0000 

0.2042 

0.2062 

0.0000 

0.0338 

0.0000 

15 

0.0000 

0.0000 

0.0000 

0.1916 

2.1228 

0.0000 

0.0154 

0.0000 

16 

0.0000 

0.0000 

0.0000 

0.0842 

0.3282 

1.0000 

0.0842 

1.0000 

Sum 

1.0000 

1.0000 

1.0000 





1.0000 

E(Xm) 

4.2848 

4.4840 

4.6854 

13.0336 

14.2047 

16.0000 

8.2938 

16.0000 

E(Yrn) 

0.1481 

0.2075 

0.2799 

8.0336 

9.2047 

11.0000 

3.5080 

11.0000 
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Table 1.54 : Probabilities P n (t) P(n) and means for m/m /r/n with r =6 X=0.08 p =0.15 
t=10, i = 0,1,2,.. ......16, and Q„ (t) with p 0 (i) =1/17 


i/n 

0 

1 

2 

14 

15 

16 

Qn ( t ) 

P(n) 

0 

0.0159 

0.0123' 

0.0095 

0.0000 

0.0000 

0.0000 

0.0038 

0.0000 

1 

0.0657 

0.0545 

0.0449 

0.0003 

0.0001 

0.0000 

0.0184 

0.0000 

2 

0.1358 

0.1199 

0.1048 

0.0015 

0.0006 

0.0000 

0.0454 

0.0000 

3 

0.1869 

0.1747 

0.1613 

0.0043 

0.0019 

0.0000 

0.0757 

0.0000 

4 

0.1920 

0.1893 

0.1839 

0.0019 

0.0042 

0.0000 

0.0958 

0.0000 

5 

0.1562 

0.1619 

0.1650 

0.0153 

0.0072 

0.0000 

0.0982 

0.0000 

6 

0.1034 

0.1123 

0.1199 

0.0212 

0.0104 

0.0000 

0.0849 

0.0000 

7 

0.0643 

0.0735 

0.0824 

0.0294 

0.0148 

0.0000 

0.0741 

0.0000 

8 

0.0377 

0.0454 

0.0537 

0.0389 

0.0201 

0.0000 

0.0653 

0.0000 

9 

0.0209 

0.0266 

0.0332 

0.0484 

0.0257 

0.0000 

0.0578 

0.0000 

10 

0.0110 

0.0147 

0.0194 

0.0561 

0.0305 

0.0000 

0.0507 

0.0000 

11 

0.0055 

0.0078 

0.0108 

0.0600 

0.0333 

0.0000 

0.0435 

0.0000 

12 

0.0026 

0.0039 

0.0057 

0.0583 

0.0329 

0.0000 

0.0356 

0.0000 

13 

0.0012 

0.0018 

0.0029 

0.0502 

0.0287 

0.0000 

0.0271 

0.0000 

14 

0.0005 

0.0008 

0.0013 

0.0364 

0.021 

0.0000 

0.0180 

0.0000 

15 

0.0002 

0.0002 

0.0006 

0.0186 

0.0108 

0.0000 

0.0088 

0.0000 

16 

0.0002 

0.0004 

0.0007 

0.5520 

0.7578 

1.0000 

0.1969 

1.0000 

Sum 


ngggs 




B ‘ 



E ( Xt ) 

4.2378 

RSI 

4.7793 


14.6366 


8.5799 


E ( Yt ) 


0.6783 

0.8278 

8.4257 

9.4257 
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Table 1.55 : Probabilities P m (n) p(n) for geom. (n) /geom.(n) /r/n for 
i -channel busy period with r =6 , N =20 , X =0.8 , p =0.1 5, m =10, i =5,6. 


i =6 


P m (n) P(n) 


n P m (n) P(n) 


0.5507 1.0000 


0.1367 0.0000 


0.2083 0.0000 


0.0798 0.0000 


0.0204 0.0000 


0.0036 0.0000 


10 0.0005 0.0000 


1 1 0.0000 0.0000 


12 0.0000 0.0000 


13 0.0000 0.0000 


14 0.0000 0.0000 


15 0.0000 0.0000 


16 0.0000 0.0000 


17 0.0000 0.0000 


18 0.0000 0.0000 


19 0.0000 0.0000 


5 0.8524 1.0000 


6 0.0108 0.0000 


7 0.0186 0.0000 


8 0.0227 0.0000 


9 0.0231 0.0000 


10 0.0208 0.0000 


11 0.0169 0.0000 


12 0.0127 0.0000 


13 0.0088 0.0000 


14 0.0057 0.0000 


15 0.0035 0.0000 


16 0.0020 0.0000 


17 0.0010 0.0000 


18 0.0005 0.0000 


19 0.0003 0.0000 


20 0.0002 0.0000 



S] 






1.6 NUMERICAL COMPUTATIONS OF DISCRETE-TIME SOLUTION 
FOR A MULTI-SERVER QUEUE WITH BALKING AND RENEGING 

We analyze a discrete time multi-server queue with balking and reneging given the 
initial state. We also discuss the case when the initial state is arbitrary. We give closed form 
solution to this class of problems in terms of roots of polynomial in z-transform and results are 
computed even when the matrices involved are large. Its is also shown, how the results in the 
continuous case can be obtained. 

This section should simulate the study of discrete time models in other areas such as 
computer science. Finally extensive numerial computations were performed in order to judge 
the accuracy of the results. Case of Machine interference problems is also given. 

Assumption : 

1 . The queue size is finite. 

2. Inter-arrival and service probabilities are dependent on the state of the system. 

3. Inter-arrival and service time distributions are geometric but independent of time. 

4. Queue discipline is first come first serve. 

Number of customer in the queue at time epoch k. 

Maximum queue size. 

Inter-arrival probability when n customer are in the system. 

e n (1 - M’n) i^nO-eJ 

Analysis: 

We develop a general discrete time Markov Model for a finite waiting space queueing 
system and analyze the effects of customer impatience on its transient behaviour. Impatience 
can be due to Balking, Reneging or both. Balking is the reluctance of a customer to join the 
queue upon arrival. Reneging is the reluctance of a customer to remain in the queue after joining 
it and leaving the queue without being serviced. The initial number of customers c will not 
renege because of their immediate entry to the service facility. Still these c customers join the 
queue with some balking probability. 

We assume tjjat inter-arrival and service times have geometric distributions with param- 
eters eand p. An arriving customer balking with probability n/N, n = 0,1, 2, 3 ....N. Thus 


Notations : 

x k •• 

N : 

e : 
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inter-arrival probability may be defined as 

e n = e (1 -n/N) ; 0 <n<N 

A customer may renege after joining the queue if hear she decides t certain waiting time 
will be larger that can be tolerated. This reneging time is assumed to have a geometric 
distribution with parameters Q. Since any one of (n - c) customers may renege, the reneging 
probability may be expressed as. 

0 for 0<n<c-l 

(n-c)Q for c < n < N 

Thus the service probability mat be expressed as 

'np for 0 < n <N 

-cp+(n-c)Q for n<c<N 

Let X k be the number of customers in the system at discrete time epoch k. Then X k , k>0 

is an integer valued discrete stochastic process taking values 0,1,2, ..,N. X k = n (0 <n<N) 

implies that there are n cutomers in the system. The process X k behaves as a discrete-time 
Markov process and represent the state of the system. 

Denote the probability that the system is in state n at the m th epoch as P m (n) (0 < n < N). 
The following difference equations may easily be written. 

P m+ , (0) 


p 
r n 


■ c 


■ m+1 
3 

" m+1 


P-.(n) 


P m+1 (c) 


P m+1 (n) 


P m (0) (l- €) + P m (l)p (l-(N-l)e/n) ....(1.33) 

= PJn)(l -(N-n)e/n-np+2nep(N-n)/N) + Pm(n-l) (1 -(n-l)p) 
(N - n + l)e/N + P m (n+1) (n+l)p (l-(N-n-l))e/N) 

1 < n < c-1 (1.34) 

= P m (c) (1 - (N-c)e/N - cp + 2cep (N-c) /N) + P m (c-1) (1 - (c-l)p) 
(N-c+l)e/N + P m (c+1) (cp+Q) {1 - (N-c-l)e/N} ...........(1.35) 

= P m (n) (1- (N-n)e/n - cp - (n-c)Q + 2 (cp + (n-c)Q) (N-n)e/N) + 
P m (n-1) (1 -cp - (n-c- 1 )Q)(N-n+ 1 ) e /N + P m (n+1) (cp + (n-c+l)D) 

{l-(N-n-l)e/N} c+l<n<N-l (1.36) 

P m (N) (1 - cp - (N-c)Q) + P m (N-l) (1 - cp - (N-l-l)Q)e/N 

with P 0 (i) - 1 , 0 < i < N ...(1.37) 

Let P z (n) be the p.g.f. of P m (n) defined as 

P z (n) = 2 z m P (n) I z | < 1 

Taking p.g.f. of equations (1 .33) to (1 .37) we get AP = [S ko ,8 kl ....... S^] 7 .,.(1.38) 

Where A is (N+l) x (N+1) ridiagonal matrix with real coefficient, p is a (N+l) x 1 


P m+1 (N) 
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column vector and 5 kl is the kronecker delta defined as 


{; 


Defining 


1/z ; k=i 

•0 ; ofterwise 

s = (1 - Z)/z we have 


A(s) 


S +T 0 -<)>, 0 

-T 0 8+^,+T, > 2 


0 

0 


0 

0 


- T N-2 S+T N-1 + Vi - ( t , N 
S+< t>N 


-T, 


0 0 0 
0 0 

Where T, and <)>, are defined above and 

P z (0) 

P z (n) 

P Z (N-1) 

P Z (N) 

From equation (1.38) using cramer rule, we may determine P z (n) explicitely as 


P z (n) = 


A „ + .( s )l 


1 A(s) | 


0 <n <N 


Where A n+i (s) is obtained from A(s) by replacing the (n+l) th column of A(s) by the 
R.H.S. in (1 .38) and | A(s) | is the determinant of A(s). 

We may observe that | A(s) | = sg N (s), Where g n (s) satisfies the recurrence relation. 
g n (s) - (S+T N x +4) N . n+ .) g n .,(s) + T N . n+J 4> N . n+1 
g.!<s) = o and g 0 (s) = 1 

g n (s) may also be expresed as the determinant & NxN real symmetric matrix g(s) as 
Hh +T 0 S i T i 0 - - - 

S+ V T . S? z- - - 


g(s) : 


I ' I 


s’/^N-Z ^N-Z S+( t ) N-l +T N-l cAn-I^N-I 

S+ V T n-l 


T 

s /T N -1 X N-1 
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The zeros of the g N (s) are the negatives of the eigen values of the matrix g(0). g(0) is a 
positive definite symmetric tri-diagonal matrix. Hence the eigen values are real, positive (>0) 

and distinct. Hence the root of g N (s) are real, negative and distinct. Let a p a 2 , be the 

roots of g N (s). Thus 


P z (n) = 


\.,(*) l 


s.n (s-a.) 


0 <n<N 


Resolving the R.H.S. of P z (n) into partial fractions, replacing s by ( 1 -z)/z and 
comparing the coeff. of z m we have. 


Case I 

For 0 < i < c-1 
1 ! 


N 


0 <n<i 


N 

P m ( n ) = b n + | J a k „( 1 - a k ) m 5 


n = i 


b n + (e/N)"- 1 n (1-jp) Z (l-a k ) m , i < n < c 


b n + "^T (e/N)"- 1 n| (1-jp) k n (l-(cp + (k - cQ)) Z ^(l-a k ) m , c<n<N 


N 


Where X n = H (1 - e + (j+l)e/N) 
Case II 


For c < i < N 

b . + j£ 5p L M.— 'ff'Cw+a-c+l-kjn). X.JJa^l-o,)-, ;c£n<c-l 


b, + Z (c|i + (i-c+l-k)Q) X n |! (!-«„)", 


; c < n < i 


P m (n) = b+Z a Kn (l-a K ) m , 


n = i 


N 


b - + -^r ( s/N )“ ' n ((i-cn+(k-c)£5) X g a* (l-oj” ; i+1 s n ^ 
Where a^’s are defined as. 


N 



With C n (s) and D n (s) being the determinants obtained by the bottom right and top left 
(nxn) square matrices farmed from A(s) such that 
|A(S)I = C N ,,(s) = D wl (s) 

C n (s) and D n (s) may be determined by the following recursive relations. 
c „( s ) = ( s+< j>M.„ + , + T_.) C ,(s) - T_, <L C Js), 1 < n < N+l 


N r (N-r)!G! p c .JI (cp+(i-c)Q) n (1-e+ie/N) 


It may be remarked that for probabilities P m (n) to remain bounded | l+a k | <1, which is 
true if <j) k + T k j < 1 . Under this condition, the sum of the absolute values of the elements in each 
row of the matrix D(0) is less than 2 and hence from the Gerschgorin’s theorem | a k | <2 (see) 
Hunter J.J. (1983), which implies |l+aj <1. 
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Using IMSL packages, we can get the eigen values (roots of g(0) and (g N (5)). 

Since (l+a K ) m -> 0 as m-»oo, the steady state distribution of P m (n) may be defined as. 
P(n)= fe p m( n ) = b n 0<n<N 

It may be noted that teh values of P m (n) have been expressed as the sum of two expres- 
sions, one pertaining to the steady state and other to the transient state. 

IMPORTANT PERFORMANCE MEASURE : 

1 . Mean number of customers in the system at epoch m 

E(X,)= f o nP» 

2. Mean number of customers in the queue at epoch m 

E (YJ = I! (n-c) P» 

3 . Probability there are r or more customer in the system at epoch m 

£ p» 

n=r 

4. Probability all servers are busy at epoch m 

E(Z ) = 2 P (n) 

5. Relaxation Time (RT) (a meassure of the length of time required by the system to settle 
to its stready state). [Morse P.M. (1958)]. may be defined as. 


N 

min (-ai) 
i=i v ' 


If m»RT then Pm(n) « p(n) 


The probability of balking at epoch m 


E(AJ=2,P)P» 

n=0 

The probability of waiting up to epoch m in the queue by those joi ning it. 


E(B m )= 2,(l-n)/NP» 


GENERAL CASE 


We have assumed that the initial queue size is fixed and equal to i i.e. 1/z occurs in only 
one position in the initial probability vector. This assumption is important when we are inter- 
ested in the transient solution. The steady state solution does not depend on the initial probabil- 


ity vector. 


When there are than one non-zero elements in the initial probability vector. The 


1 $ 


1 


lli 
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probability Q m (n) (n o,l,2 .N) defined as the probability of n customers in the system at 

epoch m irrespective of the state of the system may be defined as. 

Q m (n) = £ P m (n,i) P 0 (i) 0 < n < N 

Where P 0 (i) is the ith element of the initial probability vector and P m (n,i) is the 
probability of n sustomers in the system at epoch m assuming i as the initial number of 
customer. 

CONTINUOUS TIME CASE 

Leting e n = e n 0 + 0(0), p~ p n 0 + 0(0), m = t and m+1 = t + 0 in equations (1.33) to 
(1.37), the difference equations in m can be transformaed into differential equations in t. The 
roots equation can be changed to get the solution. The roots equation can be changed to get the 
continuous time solution from the final discrete time solution. The roots of a k of | A(s) | = 0 are 
transformed to oc k ’0. (1 + a k ’0)m tends to e k m t in continuous time when t is divided into m sub- 
interval each of length 0 such as t = m0. The parameter s itself may be treated as the transform 
parameter in continuous time. R.H.S. of (1.38) will have 1 in the ith place instead of 1/z. 

Treating Q= 0, e. = Ne, (i = 0,1,2 ,...,N), bj’s represent the steady state probabilities 

for a geom(n)/geom(n)/c/N machine interence model. 

NUMERICAL ANALYSIS 
Case (I) Balking and Reneging : 

Discrete Case : 

Assume c = 5, N = 20, e =0.5, p = 0.15, m= 10, Q = 0.005, e. = (l - i/N)e, p. = pi for 
0 < i < c, p. = pc + (i - c)Q for c < i < N and P 0 (i) = 1/21 . Table 1 .61 gives the probabilities P m (n) 
for different i(0 < i < 20), the unconditional probabilities Q m (n) and the steady state 
probabilities P(n). The last five rows give the values for E(XJ, E(Y m ), E(ZJ, E(A m ) and E(B m ). 
The epoch to reach steady state is m « 300 which is » RT = 32. 

Continuous Case 

Assume c = 5, N = 20, e = 0.5, p = 0.15 t = 10, Q = 0.005, e =(1 - i/N)e, p. = ip for 
0 < i < c p. = cp + (i - c)Q, for c < i < N and i(0 < i < 20 ) , the unconditional probabilities Q t (n) 
and the steady state probabilities P(n). The last five rows give the values for E(X t ), E(Y t ), E(A t ) 
and E(B ( ). The time to reach steady state is t « 260 which is »RT = 31. 



Case (II) Machine Interference Model 
Discrete Case : 

Assumec=5, N=20, e=0.Q4, p=0.1, m=10, e.=(N-i)e, p^ip for 0 < i < c, p^cp for 
c < i < N and P 0 (i) = 1/21 . Table 1 .63 gives the probabilities P m (n) for 0 < i < 20, Q m (n) and P(n). 
The epoch to reach steady state is m «220 which is »RT=22. 

Assume t=10 and rest of the parameters as in table 1.63, Table 1 .64 gives the probabili- 
ties P t (n) for different i(0 < i < 20), Q t (n) and P(n). The time reach to steady state is t «130 which 
is »RT=20. 

Finally the accuracy of the values (roots) that is needed in such problems is very very 
large because of the recurrence relations involved in computing the probabilties. 


Table 1.61 probabilities P m (n), Q m (n), P(n) and some important performance measures for 
Geom(n)/ Geom(n)/ c/N balking and reneging problem with C=5,N=20,m=10, e=0.5, 
p=0.05, 0=0.005, i=0,l, 2, ......... 20 and P 0 (i)= 1/21. 



0 

1 

2 

18 

19 

20 

Qm ( n ) 

P ( n ) 

0 

0.0060 

0.0025 

0.0010 ... 

0.0000 

0.0000 

0.0000 

0.0005 

0.0000 

1 

0.0480 

0.0252 

0.0124 ... 

0.0000 

0.0000 

0.0000 

0.0045 

0.0002 

2 

0.1552 

0.1019 

0.0620 ... 

0.0000 

0.0000 

0.0000 

0.0183 

0.0015 

3 

0.2673 

0.2188 

0.1643 ... 

0.0000 

0.0000 

0.0000 

0.0432 

0.0068 

4 

0.2712 

0.2780 

0.2574 ... 

0.0000 

0.0000 

0.0000 

0.0674 

0.0205 

5 

0.1673 

0.2161 

0.2468 ... 

, 0.0000 

0.0000 

0.0000 

0.0754 

0.0420 

6 

0.0654 

0.1096 

0.1580 ... 

. 0.0000 

0.0000 

0.0000 

0.0721 

0.0712 

7 

0.0167 

0.0379 

0.0711 ... 

, 0.0000 

0.0000 

0.0000 

0.0679 

0.1058 

8 

0.0026 

0.0086 

0.0219 ... 

. 0.0000 

0.0000 

0.0000 

0.0654 

0.1372 

9 

0.0003 

0.0012 

0.0045 .. 

. 0.0000 

0.0000 

0.0000 

0.0645 

0.1546 

10 

0.0000 

0.0002 

0.0006 .. 

. 0.0003 

0.0000 

0.0000 

0.0642 

0.1504 

11 

0.0000 

0.0000 

0.0000 ... 

. 0.0029 

0.0005 

0.0001 

0.0642 

0.1257 

12 

0.0000 

0.0000 

0.0000 .. 

. 0.0152 

0.0040 

0.0007 

0.0641 

0.0893 

13 

0.0000 

0.0000 

0.0000 .. 

. 0.0548 

0.0203 

0.0055 

0.0631 

0.0534 

14 

0.0000 

0.0000 

0.0000 .. 

. 0.1354 

0.0690 

0.0267 

0.0637 

0.0264 

15 

0.0000 

0.0000 

0.0000 .. 

. 0.2301 

0.1603 

0.0860 

0.0619 

0.0108 

16 

0.0000 

0.0000 

0.0000 .. 

. 0.2633 

0.2525 

0.1873 

0.05 61 

0.0034 

17 

0.0000 

0.0000 

0.0000 .. 

. 0.1939 

0.2619 

0.2725 

0.0439 

0.0008 

18 

0.0000 

0.0000 

0.0000 .. 

. 0.0843 

0.1673 

0.2533 

0.0263 

0.0001 

19 

0.0000 

0.0000 

0.0000 ... 

0.0184 

0.0570 

0.1359 

0.0104 

0.0000 

20 

0.0000 

0.0000 

0.0000 ... 

0.0014 

0.0072 

0.0320 

0.0019 

0.0000 

Total 

1.0000 

1.0000 

1.0000 .. 

. 1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

E ( Xm ) 3.6243 

4.0819 

4.5603 .. 

.15.6806 

16.4180 

17.1554 

9.9533 

9.2255 

E ( Ym ) 0.1077 

0.2166 

0.3869 .. 

. 10.6809 

11.4180 

12.1554 

5.1826 

4.2647 

E ( Zm ) 0.2523 

0.3736 

0.5029 ... 

. 1.0000 

1.0000 

1.0000 

0.8660 

0.9711 

E ( Am ) 0.1807 

0.2041 

0.2280 ... 

0.7840 

0.8209 

0.8578 

0.4977 

0.4613 

E ( Bm ) 0.0584 

0.1072 

0.1727 ... 

0.2160 

0.1791 

0.1422 

0.3339 

0.4836 
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T able 1.62: Probabilities P t (n), Q t (n), P(n) for m/m/c/n queue with t =10 and rest of the 


parameters as in table (1.61) 


Table 1.63 : Probabilities P m (n), Q m (n), P(n) and some important performance measures for 
Geom(n)/Geom(n)/C/N machine interference problems with C=5, N=20, m=10, 


€ = 

0.4, p=0.1 

i = 0,1,2. 


20, and P 0 (i) = 

1/21. 




0 

1 

2 

18 

19 

20 

Qm(n) 

i P(n) 

0 

0.0003 

0.0002 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1 

0.0065 

0.0044 

0.0029 

0.0000 

0.0000 

0.0000 

0.0008 

0.0000 

2 

0.0478 

0.0354 

0.0255 

0.0000 

0.0000 

0.0000 

0.0072 

0.0027 

3 

0.1651 

0.1346 

0.1073 

0.0000 

0.0000 

0.0000 

0.0305 

0.0161 

4 

0.2953 

0.2683 

0.2376 

0.0000 

0.0000 

0.0000 

0.0708 

0.0532 

5 

0.2789 

0.2861 

0.2837 

0.0000 

0.0000 

0.0000 

0.0956 

0.1021 

6 

0.1448 

0.1737 

0.1988 

0.0000 

0.0000 

0.0000 

0.0903 

0.1392 

7 

0.0493 

0.0727 

0.0993 

0.0000 

0.0000 

0.0000 

0.0807 

0.1624 

8 

0.0106 

0.0205 

0.0350 

0.0000 

0.0000 

0.0000 

0.0742 

0.1624 

9 

0.0013 

0.0037 

0.0084 

0.0008 

0.0001 

0.0000 

0.0715 

0.1392 

10 

0.0001 

0.0004 

0.0013 

0.0059 

0.0012 

0.0001 

0.0707 

0.1021 

11 

0.0000 

0.0000 

0.0001 

0.0269 

0.0089 

0.0019 

0.0705 

0.0638 

12 

0.0000 

0.0000 

0.0000 

0.0796 

0.0374 

0.0130 

0.0703 

0.0239 

13 

0.0000 

0.0000 

0.0000 

0.1605 

0.1022 

0.0511 

0.0694 

0.0150 

14 

0.0000 

0.0000 

0.0000 

0.2279 

0.1896 

0.1291 

0.0656 

0.0055 

15 

0.0000 

0.0000 

0.0000 

0.2293 

0.2442 

0.2191 

0.0563 

0.0019 

16 

0.0000 

0.0000 

0.0000 

0.1620 

0.2183 

0.2537 

0.0408 

0.0004 

17 

0.0000 

0.0000 

0.0000 

0.0782 

0.1328 

0.1983 

0.0229 

0.0001 

18 

0.0000 

0.0000 

0.0000 

0.0243 

0.0522 

0.1002 

0.0092 

0.0000 

19 

0.0000 

0.0000 

0.0000 

0.0043 

0.0119 

0.0296 

0.0023 

0.0000 

20 

0.0000 

0.0000 

0.0000 

0.0003 

0.0012 

0.0039 

0.0003 

0.0000 

Total 1.0000 

1 1.0000 

1 1.0000. 

.. 1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

E(Xm) 4.4846 4.7349 

5.0022 . 

..14.4807 

15.1456 

15.8104 

9.4599 

7.7358 

E(Ym) 0.2809 0.3972 

0.5432 . 

.. 9.4807 

10.1456 

10.8104 

4.6170 

2.8301 

E(Zm) 0.4851 

0.5571 

0.6266 . 

..1.0000 

1.0000 

1.0000 

0.8906 

0.9278 

E(Am) 0.2242 0.2367 

0.2501 . 

.. 0.7240 

0.7573 

0.7905 

0.4730 

0.3868 

E(Bm) 0.1406 0.1833 

0.2300,, 

.. 0.2760 

0.2427 

0.2095 

0.3654 

0.4770 




Table 1.64 : probabilities P t (n), Q t (n), (Pn) and some important performance measures for 
a m/m/c/N machine interference with t=10 and rest of the parameters as in table 1.63 



0 

1 

2 

18 

19 

20 

Qm(n) 

P(n) 

0 

0.0078 

0.0054 

0.0037... 

0.0000 

0.0000 

0.0000 

0.0011 

0.0006 

1 

0.0429 

0.0325 

0.0243 ... 

0.0000 

0.0000 

0.0000 

0.0073 

0.0044 

2 

0.1115 

0.0923 

0.0750 ... 

0.0000 

0.0000 

0.0000 

0.0228 

0.0167 

3 

0.1820 

0.1634 

0.1434 ... 

0.0000 

0.0000 

0.0000 

0.0456 

0.0402 

4 

0.2080 

0.2008 

0.1892 ... 

0.0000 

0.0000 

0.0000 

0.0657 

0.0683 

5 

0.1736 

0.1793 

0.1805 ... 

0.0002 

0.0001 

0.0000 

0.0720 

0.0874 

6 

0.1247 

0.1381 

0.1493 ... 

0.0005 

0.0002 

0.0001 

0.0746 

0.1049 

7 

0.0774 

0.0902 

0.1071 ... 

0.0015 

0.0006 

0.0003 

0.0747 

0.1175 

8 

0.0415 

0.0530 

0.0664... 

0.0040 

0.0019 

0.0008 

0.0737 

0.1222 

9 

0.0193 

0.0262 

0.0356 ... 

0.0099 

0.0050 

0.0024 

0.0725 

0.1173 

10 

0.0077 

0.0112 

0.0163 ... 

0.0225 

0.0124 

0.0063 

0.0715 

0.1032 

11 

0.0026 

0.0041 

0.0064... 

0.0457 

0.0275 

0.0153 

0.0706 

0.0826 

12 

0.0008 

0.0013 

0.0021 ... 

0.0820 

0.0546 

0.0335 

0.0693 

0.0594 

13 

0.0002 

0.0003 

0.0006 ... 

0.1281 

0.0954 

0.0650 

0.0669 

0.0380 

14 

0.0000 

0.0001 

0.0001 ... 

0.1706 

0.1441 

0.1104 

0.0624 

0.0213 

15 

0.0000 

0.0000 

0.0000 

0.1887 

0.1837 

0.1607 

0.0544 

0.0102 

16 

0.0000 

0.0000 

0.0000... 

0.1668 

0.1912 

0.1950 

0.0430 

0.0041 

17 

0.0000 

0.0000 

0.0000... 

0.1115 

0.1548 

0.1893 

0.0291 

0.0013 

18 

0.0000 

0.0000 

0.0000... 

0.0518 

0.0902 

0.1378 

0.0157 

0.0003 

19 

0.0000 

0.0000 

0.0000 ... 

0.0144 

0.0329 

0.0669 

0.0059 

0.0001 

20 

0.0000 

0.0000 

0.0000 ... 

0.0018 

0.0054 

0.0162 

0.0012 

0.0000 

Total 1.0000 1.0000 

1.0000 . 

..1.0000 

1.0000 

1.0000 

1.0000 

1.0000 


E(Xm) 4.4244 4.7178 5.0294 ... 14.5384 15.2087 15.8790 9.5082 8.0480 

E(Ym) 0.5419 0.6793 0.8460 ... 9.5384 10.2087 10.8790 4.7681 3.2672 

E(Zm) 0.4477 0.5056 0.5644 ... 0.9999 1.0000 1.0000 0.8575 0.8698 

E(Am) 0.2212 0.2359 0.2515 ... 0.7269 0.7604 0.7940 0.4754 0.4024 

E(Bm) 0.1784 0.2108 0.2456 ... 0.2729 0.2395 0.2060 0.3507 0.4235 
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CHAPTER - TWO 

DISCRETE TIME TRANSIENT SOLUTION FOR BUSY PERIOD 
ANALYSIS, DOUBLE ENDED QUEUEING AND MACHINE 
INTER-FERENCE MODEL 

2.1 DISCRETE TIME TRANSIENT SOLUTION FOR BUSY 

PERIOD ANALYSIS 

INTRODUCTION 

Busy period analysis from an integral part of any queueing system as the 
distribution of its duration is important from several point of view. It is also helpful in the 
efficient planning of the system. A busy period is initiated with the arrival of customer who 
finds the system empty and ends when the system next becomes force. However the initial 
busy period may be initiated with any number i customers in the system Busy period 
analysis concerning simple Markovian queueing system has been carried under two basic 
assumptions namely the system is taken in ‘empty’ and ‘idle’ initial condition and the 
system has infinite waiting room capacity. No attempt seems to have been made to obtained 
a transient solution for such a queueing system in discrete time. 

This chapter provides transient solution in discrete time for busy period analysis in 
Queueing theory under arbitrary initial condition and finite waiting space. It is also shown 
that the results in continuous time can be easily obtained in explicit closed form expressions 
are obtained in terms of roots. Extensive numerical results are also given using the closed 
form expression. 
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ASSUMOTION - 

1 • The queue consists of finite waiting space . 

2. Inter arrival and service propabities are independent of the system . 

3. Inter-arrival and service time distributions are geometric but independent of time . 

4. Queue discipline is first in first out (F.I.F.O) 

NOTATION- 

Xk : Number of customers in the system at epoch k 
N : Size of the queue (maximum) . 

L m : Inter-arrival probabilities when n customers are in the system. 
fi n : Service probabilities when n customers are in the system 
v|/ n :A. n (l-|i n ) (()„: finCl-^- 

ANALYSIS OF THE MODEL :- 

We define i-channel busy period 0 <i < N to begin with an arrival to the system at 
an epoch m, when there are (i-1) castomers in the system to the very next epoch (m+1), 
when there are again (i-1) castomers in the system. 

The following differential-difference equations may be defined as 

P m+1 (i-l)-Pm (i-1) = Wi ...(2.1) 

Pm+l(i) - Pm© = - (Vi+^i) Pm(i) +<j)i+l Pm(i+1) 

Pm +1 (n) - Pm(n) = - (\f n+<])n) Pm(n) +Vf/ n -l P m (n- 1 ) + (j ) n+1 P m (n+ 1 ) i<n<N- 1 

Pm + l(N) - P m (N) = - (1)N Pm + Vf/N-1 Pm(N-l) 

Let P z (n) be the probability generating function (p.g.f.) of P m (n) defined as. 

CO 

G(z,n)= Z z m P m (n); |z|<l 

m=0 


Taking the p.g.f of equation (2.1) to (2.4), for this we multiply the equation by z m and 

taking summation from 0 to co form and using (l-z)/z = s we get 

sPz(i-l) — <[>i Pz(i) ........ (2.6) 

(s+<t>i+H'i) Pz(i) - <t>i+i P Z (i+1) = 1/z (2.7) 

Pz(n-l) + (s+v|/ n +<))n) Pz(n) - (j)N+i Pz(n+1) = 1/z . .. ..... (2.8) 

-Vn-1 Pz(N-l) + (s+<j) N ) Pz(N) = 1/z ........ (2.9) 

We get 

AP = [5k(i-l), 8ki, 5k(i+i), '8 icn] / — (2.10) 

Where A is real tri-diagonal matrix of order (N-i+2) x (N-i+2) and P is a column vector of 
5ki is a Kronecker delta defined as 


Ski : 


1/z ; k=i 

0 • otherwise 

Defining s = (l-z)/z, we get 
i-1 


i-1 


i+1 


A(s)= 


N-l 

N 


l 

-Vi 


i+1 

0 


0 

0 


S+(j)i+V]/i -<j}j+l 

—VJ/j s+<j)j+l 


0 

0 


0 

0 


0 

0 


i+2 

0 

0 

-(j)i+2 


0 

0 


N-2 

0 

0 

0 


N-l 

0 

0 

0 


N 
0 
0 
0 


-Vn-2 s+<j>i n-i+Vn-i -<|>n 

0 "VN-1 S+(j)N 

(N-i+2) x (N-i+2) 
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Pz(n) = 


P*C) 

P z (i+1) 


P Z (N) | (N-i+2) x 1 
From equation (2.10), using Cramer’s rule Pz(n) are explicitly determined as 

I An-i+2 (S) | 


Pz(n) 


I A(s) | 


i - 1<N 


Where A n -i +2 is obtained from A(s) by replacing the (n-i+2) th column of A(s) by 
R.H.S. of eqn. (2.10) and | A(s) | is the determinant of A(s). 

Applying some row and column transformations on | A(s) | , it may be expressed as 
I D(s) | where D(s) is a real, symmetric, tridiagonal matrix of order 
(N-i+1) x (N-i+1). 

[s+i -A/vpi (j»i 0 ..... 0 0 0 

I -Vv)/i <J>i S+(j)i+\|/i -Vv|/ i+ i(j)i + i ..... 0 0 0 


D(s)= 


0 

0 


0 

0 


..... W\J/N-2<j>N-2 S+^n^+VN^ -^N-l^N-l 

0 S+(|)n-2+VKN-2 

(N-i+1) x. (N-i+1) 

| D(s) | are the negatives of the eigen values of the matrix D(0). It may be observed that D(0) 
is a positive definite symmetric tri-diagonal matrix and its eigen values are positive, real and 


66 


distinct. Thus the roots of the polynomial | A(s) | are real negative and distinct (one root is 
zero). 

Let at (k=o, 1 ,2, N=i+1) be the roots of | A(s) | with k=0. Then 

N-i+l 

|A(s) I = s n (s- a k ) 

k=l 

[A n .i+2(s)| 

And hence P z (n) = — ; i-l<N 

N-i+l 

s II (s-a k ) 
k=l 

Resolving right hand side of P z (n) into Partial fractions and replacing by (l-z)/z, using 
initial condition and comparing coefficient of z m , we get. 

P m (i-1) = Coff. of z m in Pz(i-l) 

I A, (s) | 

= Coff. of z m in . 

N-i+l • 
s n (s-a.j) 
j=i 

P m (i-1) = 1+ (|,i a k (i-l) 

Pm(i) = Coff. of z m in P z (i) 

|A,(s)| 

= Coff. of z m in : 

N-i+l 

s n (s - a,j) 

j =i 

Where 

P m (i) = Coff. of z m in P z (i) 

I An-2+l (s) | 

= Coff. of z m in 

N-i+l 

S n (s-aj) 
j=l 

n-l N-i+l 

= S n xRr+i Zaic„(l+a k ) m ; i <1<N-1 

r=l k=l 


Where 


C N -i ( a k) 

a k (i-l) = 

N-i+l 

«k 11 (a k -aj) 
j=l 

Cn-i (a k ) Dj (aj) 

a ki = — 

N-i+l 

Ct k n (a K - a k ) 
j=i 

j^k 

Where C„(s) and D n being the determinants obtained by the bottom right and top left 
(nxn) square matrices farmed from A(s) such that. 

I A(s) | = CN-i+ 2 (s) = Dn-h- 2 (s). 

It may be remarked that for the probabilities P m (n), i -1 < n < N to remain bounded 
I l+ct k | < 1 which is true if(j) k + t|/ k .i < 1. Under this condition, the sum of the absolute value 
of the elements in each row of the matrix D(0) is less than 2 and hence from Gerschgorin’s 
Theorem | a k | < 2; Hunter [6], which implies | l+a k | < 1. 

C n (s) and D n (s) may be determined by the following recurrence relations. 

Assuming Co(s) and D 0 (s) = 1 


Ci(s) = s+(j) N , D(s) = s 
And X = = 0 


C k (s) - (s+V|/N+l+ k + (j>N-l-k) C k -1 - lj/N+1-k (|)N+2-k C k .2 ; 2<k<N-i+2 

D k (s) = (s+\)/i + k .2 + <|> i+ic- 2 ) D k .i - y j +k . 3 (j> i +k . 2 c k _2 ; 2 <k<N-i+2 

Using the standard IMSL packages one can find the eigen values and hence the verse 
of polynomial | a(s)|. Further using recurrence relation for C k (s) and D k (s) can easily be 


evaluated. 


Using IMSL to find eigen values for large N would require much larger memory, 
hence for greater precision, for large N(N>200) one might be force to use the main frame 
computer. The first author has also the experience of working with QROOT developed at 
RMC, Canoda see [4], 

Therefore some comment on QROOT will be in order. The accuracy of the roots by 
QROOT given by | A(s) I <10 14 is not sufficient for problems under consideration because of 
the recurrence relation involved in finding the roots and the probabilities P m (n). IMSL proves 
better than QROOT for the problem considered in this section. 

Since (l+ak) m -» 0 as m-»0 the steady state distribution is given by 
P(i-l)=l; 

P(n) = 0 i < n <N 

The solution presented here is expressed as the sum of two parts one pertainting to the 
steady state and the other to the transient. 

NUMERICAL RESULTS : 

We given below the numerical results for both the discrete and continuous cases of 
the model discussed above for the sake of convenience, the results for only moderate values 
of N are given, though there were no problems even for large values of N. 

Assume r=6, N=20, k=0.8, p=0.15, m=10, i=5 & 6 

The following table gives the probabilities P m (n) and the steady state probabilities. 

• ' ' '■ . ■ . ' . • • • • \ - I 

■ . , . " . ; ' . V ^ I 

' . . ' ■ ' , ■ Q' | 

' ' ' ' . ' . ' ' ' . , ' • : ■ . ; • . , : , . ■ - "; v j 

' . ■ ' ... ' , . . . ^ ;■ ;; .• ' ■ J 

’ 

. ■ : ' '• ■' . ■ ■ , ; , . , ■' ' 

' • , ‘ • ' " ' ■ . , J 
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i =5 i = 6 


n 

Pm(n) 

P(n) 

4 

0.5507 

1.0000 

5 

0.1367 

0.0000 

6 

0.2083 

0.0000 

7 

0.0798 

0.0000 

8 

0.0204 

0.0000 

9 

0.0036 

0.0000 

10 

0.0005 

0.0000 

11 

0.0000 

0.0000 

12 

0.0000 

0.0000 

13 

0.0000 

0.0000 

14 

0.0000 

0.0000 

15 

0.0000 

0.0000 

16 

0.0000 

0.0000 

17 

0.0000 

0.0000 

18 

0.0000 

0.0000 

19 

0.0000 

0.0000 

20 

0.0000 

0.0000 


11 ■ 

Pm(n) 

P(n) 

5 

0.8524. 

1.0000 

6 

0.0108 

0.0000 

7 

0.0186 

0.0000 

8 

0.0227 

0.0000 

9 

0.0231 

0.0000 

10 

0.0208 

0.0000 

11 

0.0169 

0.0000 

12 

0.0127 

0.0000 

13 

0.0088 

0.0000 

14 

0.0057 

0.0000 

15 

0.0035 

0.0000 

16 

0.0020 

0.0000 

17 

0.0010 

0.0000 

18 

0.0005 

0.0000 

19 

0.0003 

0.0000 

20 

0.0002 

0.0000 



Finally the Discrete time models are very important for application purposes such as 
computer, it seems that this area of research has largely been ignored particularly when the 
transient solution are need. 


isi 
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2.2 DISCRETE TIME CLOSED FORM TRANSIENT SOLUTIONS FOR 
DOUBLE ENDED QUEUEING SESYTEM 

The concept of double ended queue was introduced by Kendall (1951). The 
example of such a queueing system is a taxi stand where at times passengers queue up for 
taxis and at other time taxis wait for passengers . 

Shrivastava and Kashyap (1982) obtained the transient solution of a double ended 
queue in terms of summation of the integrals of Bessel function, which is quite 
cumbersome and difficult to compute. However recently Sliarma (1990) provided a simple 
algebraic closed farm expression and easy to compute for the transient probabilities in 
continuous time seems no attempt has been made to obtain similar results in discrete time. 

It this section we obtain closed form transient solution far the double ended 
queueing system in discrete time .It is also further shown how the corresponding results in 
continuous time can be obtained. The work presented in this section is continuation with 
the earliear work done by kapur et al (1996). 

ASSUMPTION - 

(1) The numbers on the negative axis stand for taxis waiting far customers on passengers 
where as the numbers on the positive axis denote the passengers waiting for taxis 

(2) The queue consist of finite waiting space. 

(3) Probabilities of passengers arrival is k and probabilities of taxis arrival is p. 

(4) Arrival probabilities of passengers and taxis fallows geometric distributions . 

(5) The queue discipline is FCFS. 


NOTATION- 

Xk: Number of customers taxi at epoch k. 

M : Maximum number of taxis 
N : maximum number of customers . 

X : Arrival probabilities of passengers or customers, 
p : Arrival probabilities of taxi . 

<j>:X(l-p),\|/:p(l-X). 

ANALYSIS OF THE MODEL : 

Let X m be either the number of customers or taxis at discrete time epoch m then {Xk}, 
k > 0 is an integer valued discrete stochastic process. 

Taking values -M , -M+l, ..... .0,1,2,.... N-l, N , X k = k (-M < k < N ) implies that there 

are either M taxis or N customer waiting at epoch k . On each arrival or service the process 
Xk behaves as a discrete time Markov process and represent the state of the system. 


The difference -differential equation may be written as 

P m+1 (-M) = P m (-M)(l- X) + P m [-(M-l)]p (X) .........(2.11) 

P m+ i [-(M-l)] = P ra [(M-l )] [(1- X)(l-p) +Xp] + P m [-(M-2)]p (l-X) 

+ P m (-M)X (2.12) 

P m+I (k ) = P m (k )] [ (1- X) (1-p) +Xp ] + P m [(k+l)]p (1- X) +P m (k-1 )X (1-p) 

;-(M-2) < k <N-2 .......(2.13) 

Pm + i (N-l ) = P m (N-l) [ (1- X) (1-p) +Xp ] + P m ( N )p +P m (N-2 )X (1-p) (2.14) 

Pm+i (N ) = P m (N ) (1-p) + P m (N-l) A,(N-1 ) (2.15) 

With P 0 (-M) = P 0 ( i) = 0, -(M-l) < i <N 


Let P z (n) be the probabilities generating function (p.g.f) of P m (n) defined as 

00 

P z (n)= Z z m P m (n) , | z | <1 

m= 0 

Now taking the P.g.f of the equation (2. 11) to (2.1 5) we have . 

s P z (-M ) - p (1- X) P z [-(-M-l)] = 1/z ...(2.16) 

X P z (-M ) +(s+<j) -t-xp) P z [(-M-l) ] -<j> P z [-(M-2) ] = 0 (2.17) 

(j) P z (k-1 ) +(s+<j) +v(/) P z (k) -y P z (k+1) = 0 (2.18) 

(M-2) Z k € N-z 

<t> Pz (N-2 ) +(s+<(> +v(/) P z (N-l) -p P z (N ) = 0 (2.19) 

<j) P z (N-l ) +( S+p) P z (N) =0 (2.20) . 

AP= [6o(-m),5o(-m+i) -§00, Soi, 8 on] ............(2.21) 

Where A is a real tridiagonal matrix of order (N-k -2) x (N-k -2) ; P is a column 
vector and 8ok the Kronecker delta defined as 


1/z ; 

k = -M 

0 ; 

otherwise 


Pz(-M) 

P z (-M+1) 
P z (-M+2 ) 


Pz (-1 ) 

Pz(0) 

Pz(l) 

Pz (N-2 ) 
Pz (N-l) 
Pz (N ) 


(N+M+l) x 1 





Using Cramer’s rule 


An+m+i.(s) 


where An+m+i (s) is obtained from matrix A by replacing the (N+M+l ) th column 


by right hand side of (2.21) and [ A (s) | is the determinant of A(s ) 


Observe that D(0) is a positive definite symmetric tri-diagonal matrix , therefore its 


eigen values are real, positive and distinct Kijima (1992 ). Hence the roots of the 


polynomial | A (s) | are real negative and distinct (one root being zero) 


,, N+M+l) be the roots of | A (s) | with (cto = a -m = 0 ) then 


N+M+l 

I A (s) | = s n (s-ai) 


And hence 


!±M±I 


N+M+l 
s FI (s-aO 


OCN-2 - ttN+M-2, 


aN-i - ajNj+M-i an - aN+M+i 


Resolving equation (2.22) into partial fraction and replacing s by ( 1-z )/ z and using initial 


condition and comparing coefficient of z m we have 


N+M+l N+M+l N+M+l 

n C-ai )■* + FI (aj ) A +(n ( 


N+M-l 


....(2.23) 

however, if one is interested in finding the value of P m (N) under are be arbitrary initial 


condition, one may obtain the probability P m (N) as 


M+N 

d) 

Pm(N)= n — - 

i = o (-oj +1 ) 

where 


N+M+l 


(« N ™£aj(l+ajr- 

j = l 


.(2.24) 


D i ( a j) 

aj - ■ 

N+M+l 

OjII ( CCj-Oj) 
i= l,^j 
and 

M+N N+M+l 

Pm(N ° = n 1 + Sajd+aj)" 

j = 0 (-aj +1 ) j = l 


if i = N+M+l 


.(2.25) 


where D; (s) being the determinant obtained by the top left (ixi) square matrix formed from 
A(s) such that I A(s) | = Dn+m +2 (s), A,n = Xn+m+i = 0, Dj (s) is obtained by the following 
recursive relation which can be easily obtained because of tri-diagonal matrix of | A(s) | . 
Assume D 0 (s) = 1, Di (s) = S + X, D.m = 0 and XN = 0 
D;(s) = (s +(|) + v)/) Dj_ i (s) - <j) i _2 y i-i Dj.i (s) , z < i < N + M 

To get the roots a, (i = 1,2, 3 ........... N, M+l), we have used the IMSL package. It may 

be noted that we find the eigen values of the matrix D (0) and roots are negatives of the 

eigen values. The accuracy of the roots equal has been verified against the request in 
algebra viz, sum of the negatives of the roots is equal to the to the sum of the elements on 
the principal diagonal of A(0). The routines of the package are quite efficient and produce 
result to a high degree of accuracy even when the matrix size is greater so. Since (1+a,) m 
-> 0 as m -+ oo, the steady state distribution is given by 
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M+N 


p.(N)= n 1,1 

i = 0 (-o.j+1 ) 

Continuous Case 
Letting X = X (A) + 0 (A) ; 

P = p (A) + 0 (A) 

Taking m = t and m+1 = t+D, in the. difference differential equations one can 
transform equation in t we can then proceed to get the continuous time solutions, of the 
transformed equations. Alternatively, one can change the root equation and them get the 
continuous time solution form the final discrete time solution proceeding this way the roots, 
a i of | A(s) | are transform to a-, A. It is then easy to see that (l+ojA) m tends to e; a+1 in 
continuous time, where t is divided into m such interval each of length A such that 
t = m A. 

Now treating X and p as the arrival rate of passenger and arrival rate of taxi 
respectively. One can get the transient solution of the continues model. Note that the roots 
of | A (s) | were formed in s = (l-z)/z. Parameters in continues case may then be treated as 
the transform parameter. Now the R.H.S of matrix equation will have 1 in the i th place 
instead of 1/z. Analytic and numerical results in this also be verified as given in literature. 

The analogy gives the fact that discrete time model discussed in this section are 
more general and correspondingly provide results in continues time with a real case. 


N+M+l Vi N+M±l 


Pm (N) ("1) N A,(<J>) N ' +M " 1 j n ( - (X ) f+ X-*! (cXi )'' n ( <Xj - CCj ) 

^ i-1 ) i=lL j =1 j ^ i 


N+M+l 


- 


2.3 DISCRETE TIME CLOSED FROM TRANSIENT SOLUTION FOR 
GEOM (n)/ GEOM (n) /1/N MACHINE INTERFERENCE MODEL 

Benson and Cox (1951) were the first to formulate an M/M/r machine interference problem 
incorporating the operative efficiency and machine availability and they have obtained 
steady-state solution : Later on Benson ( 1951) Maritas and Xirokostas (1977) and Bunday 
and Scraton (1980) have also investigated state -state measure of effectiveness. However, 
in most practical situations where the machine interference model is applicable such as 
textile industries, robot operated factories and large computer center the study of transient 
solution are rather unavoidable if the true working of these system to be monitored over a 
small horizon of time Recently Sharma (1990) obtained simple closed farm expressions 
for state probabilities for such a system in continuous time. In this section we obtain closed 
form transient solution for Geom(n)/Geom(n)/l/N machine interference model in discrete 
time. It is also shown how the corresponding results in the continuous time can be obtained. 
Further we discussed multi -operative system wit an arbitrary initial number of machine 
interference. Extensive numerical results are also give using the also form expressions for 
multi server interference model 

ASSUMPTION: 

(1) The queue consists of finite waiting space . 

(2) Probability of a machine breakdown is Land probability of repair a machine is p. 

(3) Probability of machine are break downs and repairs follows geometric distribution 

(4) The queue discipline is FCFS. 


NOTATION: 

Xk: Number of machine breakdowns at epoch k. 

N: Maximum number of machines . 

A: probability of machine breakdown . 
p: probability of machine repair. 

A n : (N-n)Afar all n=0;l;2; .N.s 

p n : pfor alln= 1, 2, N with p 0 = 0, A N = 0 

v)/ n : A n ( 1 - p„) with \p 0 = A 0 = NX , \p N = 0 
<|) n : Pn ( 1 - A n ) with <j)o = 0 , <j)i = p [ 1-(N-1) A ] 

ANALYSIS OF THE MODEL : ' 

Let X n be the number of machine breakdown at discrete time epoch m then { Xk : k > 0 } is 
an integer valued discrete stochastic process 

Taking values 0,1,2.. ....N . Xk = k implies that there are k machine breakdowns at 

epoch k . Here arrival refers to the repair of a machine. The process { Xk } behaves as a 
discrete time Markov - Process and represents the state of the system . 

Let P m ( n ) denote the probability that the system is in the n th state at the beginning of the 
m th epoch i.e 

Pm(n) = probability that there are n machine breakdowns at a given epoch m , m > 0 . We 
further assume that breakdowns at a given epoch m, m>0. We further assume that i . 
breakdowns initially need attention at m=0. 


(2.26) 



The difference - differential equations may be written as 


Pm+i(N) = P m (N) (1-p) + P m (N-l) X (1-p) 


Let P z (n) be the probability generating function (p.g.f) of P m (n) defined as 


P z (n)= 1 z m P m (n); |z| <1 (2.3 

m=0 

Taking p.g.f. of equation (2.26) to (2.29) for this we multiply by z m and taking 


summation from 0 to °o and (l-z)/z = s, we have 


(s+NX) P z (0) - p[l-(N-l)X.] P z (l) = 0 


[(N-n+ 1 ) A,]+( 1 -p)P z (n-1) + [s+p+(N-n) (1-p)] P z (n)-p[l-(N-n-l) X] P z (n+1) = 0 


[(N-i+l)X]+(l-p)P z (i-1) + [s+(N-i) X (1-p)] P z (i)-p[l-(N-i+l) L] P z (i+1) = 1 


Whenn = 0 


The above equations can be written in matrix notation as follows, 


Where P is a column vector given as 
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P z (0) 

p z (l) 


P*(N-1) 

P Z (N) (N+l) x 1 

and Sin is a Kroneckar delta defined as 
T l/z ; ifn=l 

= i (2.37) 

°0 ; otherwise 

•V. 

Using the following notations, 

Using \j/ n = X n (l-|n n ), vj/o = 7,0 = NX, \j/ N = 0 
4*n = ft n (l'7, n ), <j)o = 0, (j) 1 = [1-(N-1)]A, 

from equation (2.36) using Cramer’s Rule we may determine P 2 (n) explicitely as 


P z (n) = 


I A n+ i(s) 


0 < n <N 


,(2.38) 


Where A n+I (s) is obtained from A(s) replacing the (n+l)* column of A(s) by the 
right hand side of equation (2.36) and | A(s) | is the determinant of A(s). 

Applying some row and column transformation | A(s) | . It may be expressed as 
s D(s), where D(s) is a real, symmetric, tri-diagonal matrix of N x N. 

I D(s) | is a polynomial of degree in N is S. It may be noted that D(0) is a positive 
definite symmetric, tri-diagonal matrix. It is well known that its eigen values are positive, 
real and distinct, one root being zero. 

We observe that D n (s) satisfies the recurrence relation 

D n (s) - (S+^/N-l + <j)N-n+l) Dn-l(s) + l|/N-n+l <(>N-n+l D n -2(s) = 0 ; 1 < n < N 




1 


l 
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Let ai, 0 C 2 , cxn be the roots of D n (s). Thus 

|A n +i(s) | 

P z (n) = ; 0 <n <N (2.39) 

N 

FI (s - aO 

i=l 

Resolving the right hand side of P z (n) into partial fractions, by (1 -z)/z and 
comparing the coeff. of z m , we have 


i N 

Pm(n) =b n +l/n!p- n riiiB n Ea mn (l-a k ) m ; n=0 (2.40) 

k=l k=l 

i n 

Pm(h) — b n ■+■ n M' B n ^kn (l-a k ) 1 < n < i (2.41) 

k=l k=l 


i-1 

Where B„ =11 {l-(N-j-l)A, 
j=n 

n 

P m (n) =b„ + n a kn (l-a k ) m ; n = i ..(2.42) 

k=l 

n-1 N 

P m (n) =b„+{(n-i)!/(N-n)!}X n - 1 n (1-p) I a kn (l-a k ) m ; i+1 <n<N ...... .(2.43) 

k=l k=l 

Where a kn ’s are defined as 

C N -i (oc k ) D n (a k ) 

a kn = — , 0 < n < i 

N 

oc k II (a k -aj) 
j=l*k 

C N . n (a k ) Dj (a k ) 

a kn = , 0 < n < N 

N 

a k n (a k -ctj) 

j=l*k 
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With C n (s) and D n (s) being the determinants obtained by the bottom right and top 
left (nxn) square matrices formed from A(s) such that 
| A(s) | = Cn+i(s) = Dn+i(s) 

C n (s) and D n (s) may be determined by the following relation. 

C n (s) = (s+<j>N-n+l + H^N-n+l) C n -i(s) — Ij/N-n+l <j>N-n+2 C n - 2 (s) 5 1 ^ n < N+l ......(2.44) 

With Co(s) = 1, C.i(s) = 0 

D„(s) = (s+<j>n-i + V|/n-i) Dn-i(s) - <t> n -i ¥n-2 D„. 2 (s) ;l<n<N+l (2.45) 

With Do(s) = 1 and D.,(s) = 0 
C ' r ' 




n n {i-(-i)p} 

l+Z N C r (X/p) r i=l 

r=l 

r 

n {i-(N-oM 

i=i 


n=l 


V 


1 


1 + N (X/p) 


N 

{1-(N-1)X}S 


N! Ar (1 - jj.) 1 


r-1 


r 


r=2 r 

N r (N-r)! (p) r ri{l-(N-i)X} 
i=l 


n {1 -(-1) |x} 

N C r (>7p) r i=l 


i <r< 1 


brb 0 ' 


i 


n {1 -(N-i) X} 

i=l 


N!/(N-r)!(X/Np) f 


0 -lt) r 


n{i-(N-i)A> 

i=i 


2 <r <N 


.(2.46) 


.(2.47) 



It may be noted that for probabilities P m (n) to remain bounded if 1 1+a.k I <1 which 


is true if (<j>k + \|/ k .i)<l, under this condition the sum of the absolute value of the elements in 


each row of the matrix D(0) < 2. 

Hence | ak I < 2 = 1 1+ock | <1 


Using IMSL package, we can get the eigen value of D(0). The routiences of this 


packages are quite efficient and produce result to a high degree of accuracy even when the 


matrix size is large that is > 50 


Since (l+ak) m ->0 as m->oo the steady distribution Pm(n) may be defined as 


Discrete Transient Solution For Geom(N)/Geom(N)/C/N Multiserver Machine 


Interference Model 


This model is an extension of the discrete transient solution of machine interference 


problem discussed above, here instead of single server there are c(l < c < N). Under 
indentical assumption made for a single server, we have for this model. 


We assume that inter-arrival and service times have geometric with parameters X 


and p respectively. 
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Denoting the probability that the system is in state n at the m* epoch as P m (n) 

(0 < n<N). 

The following difference equations may be easily written as. 

P m+ i(0) = P m (0) (1-NA.) + P m (l)p [ 1 -(N- 1 ) X] ....... (2.48) 

Pm + i(l) = Pm(l) (1-p) [ 1 -(N- 1 ) X] + P m (0) m + P m (2)p [l-(N-2j X] (2.49) 

Pm+i(n) = P m (n) [l-(N-n) X] (1-np) + P M (n-l){(N-n+l)X} {l-(n-l)p} + P m (n+1) 

[l-(N-n-l) X] (n+1) (a ; 1 <n<c ..........(2.50) 

P m+ i(c) = P m (c) [l-(N-c) X] (1-cp) + P m (c-1) [l-(c-l) p] (N-c+1) X + P m (c+l)(cp) 

[l-(N-c-l) X] (2.51) 

P m+ i(n) = P m (n) [l-(N-n) X] (1-cp) + P m (n-1) (1-n+l) X (1-cp) + P m (n+1) (cp) 

[l-(N-n-l) X] ; c<n<N-l (2.52) 

Pm+i(N) = P m (N) (1-cp) + P m (N-l)X (1-cp); (2-53) 

Taking p.g.f. of equation (2.48) to (2.53), for this we multiply by z m and summing 
from 0 to oo, after simplification, we have 

(s+NA.) P z (0) - p[l-(N-l) A.] P z (l) = 0 (2.54) 

-m P z (0) + [s+p+(N-l) X (1-p)] P z (l)- p[l-(N-2) X] Pz(2) = 0 ....... (2.55) 

-(N-n+l)A, [i-(n-l)] P z (n-1) - [s+np+ (N-l) X (1-np)] P 2 (n) - (n+1) p[l-(N-n+l) a] 

P z (n+1) = 0: 1 < n < c (2.56) 

-(N-c+1) A. [l-(c-l) p] Pz(c-l) + [s+cp+(N-c) X (1-cp)] 

Pz(c) -Pz(c+l)cp[l-(N-c-l) X] = 0 ....... (2.57) 

-(N-n+1) X (1-cp) P z (n-1) + [s+cp+(N-n) X (1-cp)] 

P z (n)-cp[l-(N-n-l) X] P 2 (n) (n+1) = 0 ..(2.58) 
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-X (1-cp) P Z (N-1) + (s+cp) P Z (N) = 0 
Using the notations 

(N-n)X (l-n[i) 


(2.59) 


tyn A, n (l-|in) “ a 


*|*n Pn (l“X n ) 


if 0 < n <C 

(N-n)X(l-qx) ifn>c 

With \j/ 0 = X„ = NX and vj/ N = 0 Since X N = 0 

np [l-(N-n) X] if 0 < n <c 

cp [l-(N-n) X] if c < n <N 

<|>n = Pn (I-Xn) - c|i Since Xn = 0 
(s+h/o) Po(0) = v|/i P z (l) = 0 
-Wo P z (0) + (s+p + ¥i) Pz(l) = 0 

-W„-1 Pz(n-1) + (s+v|/„+Pn) - 4*2 Pz(2) Pz(n) - 4> n +i Pz(n) = 0 
-VPc-l Pz(c-l) + (s+Cp+VJ/c) P Z (C) - 4>c+l P Z (C+1) = 0 
-H/n -1 Pz(n- 1 ) + (s+cn+v|/ n -c|>c+i) Pz(n+ 1 ) = 0 
-Wn-i Pz(N- 1 ) + (s+cp) P Z (N) = 0 
WithP 0 (i)=l; 0<i<N 

The above equations can be written, in matrix notation as 

AP = [Sjo, 5ii, M' 

Where P is a column matrix of order (N+l) x 1 
Pz(0) 


(2.60) 


(2.61) 


(2.62) 


(2.63) 


2.64 


(2.65) 


( 2 . 66 ) 


2.67) 


( 2 . 68 ) 


PzO 


(2.6 


P Z (N-1) 

Pz(N) 


(N+l) x 1 




From equation (2.68), using Cramer’s Rule, we determine P z (n) explicitely as 
I A n+ i(s) | 

P 2 (n) = , 0 < n < N ..(2.70) 

. |A(s): I 

Where A n +i(s) is obtained form A(s) by replacing the (n+l) 111 column A(s) by the 
R.H.S. of equation (2.68) and | A(s) | is determinant of A(s). 

We may observed that |A(s) | — s D(s), where D(s) satisfies the following 
recurrence relation. 

D n (s) - (s+ \j/N-n+l +<|>N-n+l) D n .i(s) + VJ/N-n+1 <t>N-n+I D n - 2 (s) = 0; 1 < n < N (2.71) 

With D-i(s) = 0 and Do(s) =1 

The zero’s of Dn(s) are the negatives of the eigen values of the matrix D(0), D(0) is 
a positive definite symmetric tri-diagonal matrix. Hence the roots of D N (s) are real negative 
and distinct. 

Let oci , a 2 ,. , (Xn be the roots of D N (s). Thus 

|A n+1 (s) | 

P z (n) = ; 0 <n <N (2.72) 

N 

s II (s - aO 
i=l 

Resolving the R.H.S. of P z (n) into partial fraction replacing s by (l-z)/z and 
comparing the coefficient of z m , we have 

N 

P m (n) =b„+(i!/n!)(p) 1 - n B n E a kn (l-a k ) m ; 0<n<i ..(2.73) 

k=l 


P m (n) = bn + Z a kn (l-a k ) m ; n = 1 (2.74) 

k=l 
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n-1 n-CN 

Pm(n) = b n + { (n-i) ! / (N-n) ! } (X/N)"' 1 U (1-jji) E a kn (l-a k ) ra ; i <n<c ....... (2.75) 

j=l k=l 


c-1 


n-c N 


Pm(n) = b n + {(n-i)!/ (N-n) ! } (A7N) n " U (l-j|i)(l-c^) I a kn (l-a k ) m ; c<n<N 


J=i 


k=l 


.(2.76) 


i-1 

Where B„ = n {l-(N-j-l)X} 
j=n 

For c < i < N 


N 


Pm(n) = b„ + { (c- 1 )!/n! } (j-i) c_n-1 (cpf* 1 B„ E a kn (l-a k ) n 

k=l 


Pm(n) 


N 

= bn + (cy)" 1 Bn E aicn (l-ct k ) m ; c<i 
k=l 


N 


P m (n) =b„ + Ea kn (l-a k ) m ; 
k=l 


n = i 


(2.77) 


.(2.78) 


.(2.79) 


m 


I 

j 


N 

P m (n) = b n + {(n-i)!/(N-n)! j (X7N) n " E a kn (l-a k ) n 

k=l 

Where a kn ’s are defined as 

Cn- 1 (otk) Dn ( a k) 


&kn 


akn 


N 

a k n (ct k - ocj) 
j=l.*k 

Cn-i (ctk) Di ( a k) 
N 

a k n (ctk - ctj) 

j=l=k 


i+1 < n < N .......(2.80) 




if 0 < n < i 


y ..(2.81) 


if 0 < n < i 


i . 1 


I 


m 
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With C n (s) and D n (s) being the determents obtained by the bottom right and top left 

square matrix formed from A(s) such that 

I A(s) | = Cn+i(s) = D n+ i(s) 

C n (s) and Dn(s) may be determine by the following recurrence relations. 

C n (s) = (s+ 4»N-n+i +M^N-n+i) C„.i(s) + vj/ N .„+i <j) N . n+ i C n . 2 (s) = 0; 1 < n < N+l .. (2.82) 

Dn(s) = (s+ (j) n -i +v|/„-i) Dn-i(s) - i|/ n . 2 <)>„., D n . 2 (s) ; 1 < n < N+l ....... (2.83) 

With Do(s) = 1 and D.i(s) = 0 


bo ’ 1 = 


N n { 1 — (i- 1 ) pj 

1 + £ N C r (X/Np) r i=l 

r=l r 

II {1 - (N-i) A,} 
i=l 


c II { 1 — (i- 1 ) c N! 

1 + S N C r (X/Np) r M + S (A/N) r 

r=l r i=c+l (N-r)!c! 

n {1 -(N-i) 
i=l 


np-Mii) 

io-cn/cn)'^ (l/(n Q i=l 

r 

n {1 — (N-i)A,} 

i=l 


.(2.84) 


P 


! 



It may be noted that the probabilities P m (n) to remain bounded 1 1+ctk | < 1 which is 
true if (\j/ k -i +<j>k) < 1, under this condition the sum of the absolute values of the elements of 

< 2. Hunder (1983) which 


the matrix D(0) < 2. Hence from Gerscehgorin’s Theorem 


implies that 1 1+ak I < 1. 

Using IMSL package we can get eigen value of D(0). Since (1+ock ) ra as m— »oo the 
steady state distribution P m (n) may be defined as 

PfnY— lim P-Jn! = b n ; 0<n<N 


Further we observe that the values of P m (n) have been expressed as sum of two 


expressions namely steady state and transient. 


GENERAL CASE: 

So for we have assumed that initial queue size is fixed and equal to i that is 1/z 
occurs only one position of the initial probabilities vector. This assumption is important 
when we are interested in the transient solution. The steady state solution does not depend 


on the initial probabilities vector. 
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We now consider a more general case of this problem, when there are more than one 
non zero element in the initial probabilities vector. The probabilities Q m (n), {n=0,l,....N} 
defined as the probability of n customers in the system at epoch m irrespective of the state 
of system may be defined as 

N 

Qm(n) = IPm(n,l) P 0 (i); 0 < n <N --(2-86) 

i=i 

Where Po(i) is the i th element of the initial probability vector and P m (n,i) is the 
probability of the n machine in the system at epoch m assuming i as the initial number of 
machine. 

CONTINUOUS TIME CASE OF GEOM(N)/C/N DISCRETE MODEL : 

Letting . = X. n 0 + 0(0) 

(J-n = p n 0 + 0 (0), 

and m = t, m+1 = t = 0 

Substitute in equation (2.48) to (2.53), the difference equation in m can be transform to 
differential equations, in t, we can proceed as above to get continuous time solution from 
the transformed equations. 

Alternatively the root of the equation can be changed to get the continous time 
solution from the final discrete time solution, the roots a k of the | A(s) I = 0 are transform 
into. It is easy to see that ( 1 +ajce) m —> subintervals each of length such as t — m0. 

Moreover parameter s may be treated as transform parameter in continuous time, the 
right hand side of eq. (2.58) will have one in the i th place of instead of 1/z. 



Important Performance Measures of Geom(N)/Geom(N)/C/N machine interference 


Discrete Model 


Using explicit expressions for P m (n), some important measures can be defined as 


under (for fixedi) 


1 , Mean number of customer in the system at epoch m 


Mean number of customer in the queue(excluding those in service) at epoch m 


Probabilities that there are r or more customer in the system at epoch m is 


4. Probability all servers are busy at epoch m 


Relaxation Time (RT) [a measure of the length of time required by the system to 

settle to its stady-state Morese (1958)] may be defined as. 


min (-aO 


if m>RT then P m (n) ~ P(n) 


6. The Probability of balking at epoch m 
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7. The Probability of waiting upto epoch m in the queue by those joining is 
N 

E(B m ) = £ (1-n/N) P m (n) 

n=c+l 

NUMERICAL RESULTS : . 

We give below the numerical results for both the discrete and continuous cases for 
the model discussed above 
Discrete Case : 

Assume c = 5, N = 20, X= 0,04, p= 0.1, m = 10 
X = (N-i), p-i = ip for 0 < i < c 
p,i = c(i; c < i < N 
and P 0 (i) = 1/21 

Table 1 gives the probabilities P m (n) for 0 < i < 20, the unconditional probabilities 
Q m (n) and the steady-state probabilities P(n). The last five rows give the values for 
E(X m ), E(Y m ), E(A m ) and E(B m ). 

The epoch to reach stady-state is m 220 is > RT = 22. 

Continuous Case : 

Assume t =10 and rest of the parameters as in Table 1 and Table 2 gives the 
probabilities Pi(n) for different i(0 < i < 20) the unconditional probabilities Qi(n) and the 

steady-state probabilities P(n). 

The time to reach steady-state is t = 130 which is » RT = 20. 

Finally the accuracy for the eigen values or roots that needed in such types of 
problems is very very large because the recurrence relations involved in computing the 

probabilities. 
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Table 1. : Probabilities P m (n), Q m (n) and P(n) and important performance measure tor 
symmetric and Geom(n)/C/N machine interference problem with C=5, M=30 and p = 0.1. 
P 0 (i) = 1/2!; i = 0,1,-.. ......... .....20 


0 

1 

2 

0.0005 

0.0005 

0.0005 

0.0059 

0.0056 

0.0053 

0.0328 

0.0314 

0.299 

0.0990 

0.0954 

0.0915 

0.1689 

0.1642 

0.1589 

0.1939 

0.1909 

0.1874 

0.1828 

0.1827 

0.1824 

0.1424 

0.1450 

0.1477 

0.0920 

0.0957 

0.0997 

0.0492 

0.0524 

0.0561 

0.0217 

0.0238 

0.0262 

0.0079 

0.0089 

0.0102 

0.0024 

0.0027 

0.0032 

0,0005 

0.0007 1 

0.0008 

0.0001 

0.0001 

0.0001 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.000 

1.000 

1.000 

6.6264 

6.6931 

6.7685 

1.8618 

1,8702 

1.9377 

0.8618 

0.8671 

0.8728 

0.3313 

0.3347 

0.3384 

0.4291 

0.4337 

0.4385 


18 

19 

20 

Qm ( n ) 

P(n) 

0.0000 

0.0000 

0.0000 

0.0002 

0.0002 

0.0001 

0.0000 

0.0000 

0.0022 

0.0027 

0.0006 

0.0004 

0.0002 

0.0130 

0.0161 

0.0031 

0.0020 

0.0013 

0.0423 

0.0532 

0.0099 

0.0069 

0.0047 

0.0801 

0.1021 

0.245 

0.0184 

0.0135 

0.1088 

0.1392 

0.0516 

0.0415 

0.0326 

0.1299 

0.1624 

0.0910 

0.0779 

0.0653 

0.1387 

0.624 

0.1342 

0.1221 

0.1089 

0.1341 

0,1392 

0.1657 

0.1597 

0.1512 

0.1176 

0.1021 

0.1706 

0.1740 

0.1746 

0.0932 

0.0638 

0.1461 

0.1573 

0.01669 

0.0656 

0.0339 

0.1031 

0.1170 

0.1311 

0.0403 

0.0150 

0.593 

0.0709 

0.0837 

0.0209 

0.0055 

0.272 

0.0344 

0.0427 

0.0090 

0.0017 


0.0097 0,0130 0.0170 0.0031 0.0004 

0.0027 0.0037 0.0051 0.0008 0.0001 

0.0001 0.0001 0.0001 0.0001 0.0001 

0.0001 0.0001 0.0001 0.0000 0.0000 

0.0000 0,0000 0.0000 0.0000 0.0000 

1.000 1.000 1.000 1.000 1.000 

10.5867 10.8801 11.1736 8.4978 7.7358 

5.5912 5.8829 6.1753 3.5737 2.8301 

09963 0.9976 0.9985 0.9422 0.9278 

0.5293 0.5540 0.5587 0.4249 0.3868 

0.4603 0.4489 0.4366 0.4679 0.4778 


Iffjfc 





j . ■ 

I 


i i ':■■■ 

i 


'% h 

"ii 


: •* 

. * '• 




Table 2. : Probabilities P t (n), Q(n) and P(n) and important performance measure for 
symmetric and M/M/C/N machine interference problem with t=10, N=20, m = 10 and 


1 = 00.04, p = 0.1, P 0 (i) = 1/2!; i = 0,1, ...20 



0 

1 

2 


18 

19 

20 

Qm(n) 

P(n) 

0 

0.0078 

0.0054 

0.0037 

... 

0.0000 

0.0000 

0.0000 

0.0011 

0.0006 

1 

0.0429 

0.0325 

0.0243 

... 

0.0000 

0.0000 

0.0000 

0.0073 

0.0044 

2 

0.1115 

0.0923 

0.0750 

* • * 

0.0000 

0.0000 

0.0000 

0.0228 

0.0167 

3 

0.1820 

0.01634 

0.1434 

... 

0.0000 

0.0000 

0.0000 

0.0456 

0.0402 

4 

0.2080 

0.2008 

0.1892 

... 

0.0000 

0.0000 

0.0000 

0.0657 

0.0683 

5 

0.1736 

0.1793 

0.1805 

... 

0.0002 

0.0001 

0.0000 

0.0720 

0.0874 

6 

0.1247 

0.1381 

0.1493 

• •• 

0.0005 

0.0002 

0.0001 

0.746 

0.1049 

7 

0.0774 

0.0920 

0.1071 

... 

0.0015 

0.0006 

0.0003 

0.0747 

0.1175 

8 

0.0415 

0.0530 

0.0664 

... 

0.0040 

0.0019 

0.0008 

0.0737 

0.1222 

9 

0.0193 

0.0262 

0.0356 

... 

0.0099 

0.0050 

0.0024 

0.0725 

0.1173 

10 

0.0077 

0.0112 

0.0163 

• • • 

0.0225 

0.0124 

0.0063 

0.0715 

0.1032 

11 

0.0026 

0.0041 

0.0064 

... ' 

0.0457 

0.0275 

0.0153 

0.0706 

0.0826 

12 

0.0008 

0.0013 

0.0021 


0.0820 

0.0546 

0.0335 

0.0693 

0.0594 

13 

0.0002 

0.0003 

0.0006 


0.1281 

0.0954 

0.0650 

0.0669 

0.0380 

14 

0.0000 

0.0000 

0.0001 


0.1706 

0.1441 

0.1104 

0.0624 

0.0130 

15 

0.0000 

0.0000 

0.0000 


0.1887 

0.1837 

0.1607 

0.0544 

0.0102 

16 

0.0000 

0.0000 

0.0000 

... 

0.1668 

0.1912 

0.1950 

0.0430 

0.0041 

17 

0.0000 

0.0000 

0.0000 

... . 

0.1115 

0.1548 

0.1893 

0.0291 

0.0013 

18 

0.0000 

0.0000 

0.0000 

... 

0.0518 

0.0902 

0.1378 

0.0157 

0.0003 

19 

0.0000 

0.0000 

0.0000 

. . . 

0.0144 

0.0329 

0.0669 

0.0059 

0.0001 

20 

0.0000 

0.0000 

0.0000 

. . . 

0.0018 

0.0054 

0.0162 

0.0012 

0.0000 

Total 

1.0000 

1.0000 

10.000 

. .. 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

E(X m ) 

4.4244 

4.7178 

5.0294 


14.5384 

15.2087 

15.879C 

9.5082 

8.0480 

E(Y m ) 

0.5419 

0.6793 

0.8460 


9.5384 

10.2087 

15.789C 

9.5082 

0.0480 

E(Z m ) 

0.4477 

0.5056 

0.5644 

. . .. • ' ’ 

0.9999 

1.0000 

1.0000 

0.8575 

0.8698 

E(A m ) 

0.2212 

0.2359 

0.2515 


0.7369 

0.7604 

0.7940 

0.4754 

0.4024 

E(B m ) 

0.1784 

0.2108 

0.2456 


0.2729 

0.2395 

0.2060 

0.3507 

0.4235 
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CHAPTER - THREE 

DISCRETE TIME TRANSIENT ANALYSIS OF A SYSTEM WITH 
QUEUE DEPENDENT SERVERS 

This chapter of thesis deals with diserete time analysis of an infinite capacity 
Markovian Queueing model with number of heterogeneous channels depending on system 
size. The model consists of one regular channel and a maximum of C-additional channel by 
defining probability generating functions. 

Introduction : 

Most of the Queueing theary literature concentrates on finding the steady state 
solutions or approximations. Very little seems to have been done to evaluate the transient 
solutions. Even at time steady state solutions are difficult to compute. 

Chaudhry, Agarwal & Templetion (1992) have mostly concentrated on this using 
the techniques of roots. Earlier attempts at finding transient solutions can be attributed to 
Tacka’s (1962) & Morse (1959). However there are computations difficulties with their 
methods. Now with the increased skill available in computations with the use of computers 
researches especially in computer science have started looking for transient solutions and 
easy to compute closed form solutions. Sharma & Dass (1988) have provided transient 
solutions to a class of Markovian Queueing models in Queueing Theory. However they did 
not concentrate on the computational difficulties of finding roots or eigen values if the 
matrices involved are large. 

Most application of Queueing theory involves queues which are emptied and 
restarted periodically. These queues never reach equilibrium state. Hence steady state 
solutions are not always adequate and it is desirable to have time dependent, solutions. 


| 

i 

i 
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Sharma and Shobha (1988) have considered multiserver Markovian Queueing model with 
finite waiting space and derived transient results. 

Sharma and Tarabia (2000) have obtained the transient state probabilities of a Singh 
server Markovian Queue with finite source in a closed form using Laplace transform 
technique. Tarabia (2001) has analysed the transient behaviour of a non-empty M/M/l/N 
queueies model using Laplace transform technique. 

In real situation it is common that the executive concerned is obliged to offer special 
service facility to his customers to gain their good will. Systems with additional servers 
may be required in places like banks, tax offices, transport facilities etc. hence a system 
with queue depended servers is proposed, Gross and Harries (1985) have mentioned that the 
transient solution of single channel queue becomes more complicated with the restriction on 
waiting room capacity is relaxed. However in this chapter we presents the discrete time 
transient solution of an finite capacity queuing model with additional servers. Here we 
assume that the system capacity is N. 

Moreover, no attempt seems to have been made to obtain similar results in discrete 
time for finite waiting space problems in Queueing theoiy. 

As the transient solution is not independent of the initial state of the system’s 
behaviors, further, some systems may not exist long enough, to reach their steady state. 

There are several systems which operate at discrete times see Kobayashi (1983). As 
a result it becomes important to study them, in such cases events are clock controlled. 

In this chapter we analyze a discrete time transient analysis of a Markovian 
Queueing model of a system with queue dependent servers. Such problems occur not only 
in Queueing theory but also in bio- science and now in computer science we given closed 
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form solution to this class of problems in terms of the roots of a polynomial in Z- transform 
and results are computed even where the matrices involved are large. It is also shown how 
the results for the continuous can be obtained. Interesting analogy exists between the 
discrete time models and their continuous time counter parts. 

Results presented in this chapter further unify the treatment given by Chaudhry, 
Kapur, Templeton (1991). It is worth noting though, continuous time models are particular 
cases of discrete time models, yet this area of research has remain neglected. Assume that 
the customers arrive individually at a service facility in accordance with a Poisson process 
having rate X. The queue discipline is FIFO. The service facility consists of one regular 
service channels. The system starts with a regular service channels. The system starts with a 
regular service channel and c- additional service channel having service rate p. 

The regular channel is always open irrespective of queue length. If the system is 
empty then the regular channel is idle. If the number of customers in the system increases to 
r (0) =1 then the regular channel will become busy. If the number of customer in the system 
is r (1) then one additional service channel with service rate p (1) is started, which will be 
dropped at the termination of a service, if the system size becomes less than r (1) at that 
epoch. 

When two channels are operating, if the number of units in the system increases 
to r (2) ( > r(l)), the second additional channels is limited to a maximum of c and when all 
the c- channels are operating together with the regular channels. 

The queue is allowed to grow without limit initially system is assumed to have 
c- units. Let p m (n) be the probability of n customers in the system at m th epoch or time slot, 
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The following difference equation may easily be written as 

p m+1 (0) -pm(0) = - Xp m (0) + pp m (l) (3.1) 

P m+i (n) - p m (n) = -(?v+(a k ) p m (n) + p m (n-1) X+ p k p m (n+1) . . ...... (3.2) 

; r (k) < n < r (k+1) 


P m +1 (r(j) -1) - Pm (r 0-1 ) = - 0+Pj-l) Pm(r(j)-1) +Pj Prn (r (j)) Pm (r (j) -2) 


Where j = 1,2, 3 (3.3) 

p m+i (n) - Pm (n) = - X p m (n-l) +p c p m (n+1) — 0+|i c ) p m (n) (3.4) 


;n>r (c) 

Defining the probability generating function as 
r (c) oo oo oo 

P(z,m,n) = S Z p m (n) z m + Z Z p m (n+r (c)) z m | z | < 1 

n = 0 m=0 m=0 n=l 

oo r(c) oo oo oo 

= Zz m z Pm(n)Z z m R(m) + z Z m Zp m (n+r-l) 
m = 0 m=0 m=0 m=0 n=l 

using the initial condition, we have 
P(z,0) =Po(n) = z K(i) orz 5(i) 

now taking the probability generating function of equation (3.1), (3.2), (3.3), (3.4) 



r(c) 

R(m) 

= Z Pm (n) 


n=0 


00 

P(z,m.n) 

= R(m) + Z pm (n+r(c)) z n 


n =1 


oo r(c) oo oo 

P(z,m,n) 

= Z Z Pm (n) + z Z Pm (n+r (c)) z n 


m=0 n=0 m=0 n=l 


P(z,m,n) = £ £ z m p m (n) 

m=0 n=l 
N 

= £ z° p 0 (n) + z 1 pi (n) + z 2 p 2 (n) +. . .. 

n=l 

Pm (n) = Coeff of z m in the expansion of P (z, m, n). 

Multiplying eq. (3.1) by z m and summing from 0 to oo, We have 

z m Pm+i (0) - p m (0) z m = -X Pm (10) z m + p Pm (1). z m 

00 00 oo 

£ Z m p m+1 (0) - p m £ (0) z m = - X Pm (0) I z m + p p m (l) . z m 
m=0 m=0 m=0 

00 

(1/z) £ z m+l p m +i (0) - p(z, m, o) = -X p(z, m, o) + p p(z, m, 1) 
m=Q 

(1/z) [p(z, m, o) - p 0 (0)] - P(z, m, 0) = - X P (z, m, o) + p P(z, m, 1) 

(1 - z)/z P (z, m, o) - 1/2 p 0 (0) = - X P(z, m, o) + p P(z, m, 1) 

Denoting (1 - z)/z = s, we have 
s P(z, m, o) - 1/z p 0 (O) = - X P(z, m, o) + p P(z, m, 1) 

(s+ X) P (z, m, o) - p P(z, m, 1) = 1/z p 0 (0) ............(3.5) 

Multiplying eqn (3.2) by z m and summing over 0 to oo, we have 
00 00 00 00 

£ p m+ i (n)z m - £ p m (n) z m = - (A,+p k ) £ P m (n) z n + 1 £ p m (n-1) z m + p R p m (n+1) z ra 
m=0 m=0 m=0 m-0 

00 00 00 

( 1/z) £ p m+1 (n) z m = 1 -£ p m (n) z m = - (X+p k ) £ p m (n) z m + X £ p m (n-1) z m + p k 
m=0 m=0 m=0 

00 

> £ Pm(n+1) z m 

m=0 

(1/z) [P(z, m, o) - P 0 (n)] - P(z, m, n) = - (X+p k ) P(z, m, n) + P(z, m, n-1) + p K P(z, m, n+1) 


— J 

. ■' '. " ' ' 

; -| 

(l-z)/z P(z, m, n) - 1/z po(n) = - 0+p k ) P(z, m, n) + X P(z, m, n-1) + p k P(z, m, n+1) 
sP(z,m,n) + P(z, m, n) - X P(z, m, n-1) - p K P(z, m, n+1) = 1/z Po(n) 

,* ■ . (s + X + p r ) P (z, m, n) - X P (z, m, n-1) - p,.P(z, m, n+1) = 1/z Po(n) ; r(r) < n < 

r( r+ 1 ) ....(3.6) 

Multiplying eq n (3.3) by z m and summing from 0 to oo, we have , 

00 00 00 00 

Z Pm +1 (r© - 1) z m - Z Pm (r(j) - 1) Z m = - (X +Pj) Zp(r(j)- l)z m + pjZpm(r(j))z m 
m=0 m=0 m=0 m=0 

°° 

+ ^Zpm (r0 - 2) z m 
m=0 

00 00 00 

1/z Z Pm+i(r(j) - l)z m+1 - Z Pm (r(j) - 1) Z ra = - (X + (Ij.i) Z Pm(rG) -1) z m + 
m=0 m=0 m=0 

GO 00 

PjZpm(r(j))z m + 1 + Z pm (r(j) - 2) z m 
m=0 m=0 

1/z [P(z, m, (r(j) - 1)) P 0 (r(j) - 1)] — P (z, m, r(j) - 1) 

= - (X + pj.i) P(z, m, r(j) - 1) + pj P(z, m, r(j)) + X P(z, m r(j) - 2) 

(1 - z)/z P (z, m, r(j) - 1) - 1/z po (r(j) - 1) 

= - (X + pj-i) P(z, m, r(j) - 1) + pj P(z, m, rQ)) + X P(z, m r(j) - 2) 
s P(z, m, 0) - 1) + - (X + pij.i) P(z, m, r(j) - l)-pj P(z, m, r(j)) - X P(z, m r(j) - 2) 

= l/zp 0 (r(j)-l) 

(s+ X + pj) P(z, m; r(j) - 1) + Mj P(z, m, r(j)) + X P(z, m, r(j) - 2) = 1/z p o (r0-l). • • • • • • -(3.7) 

Multiplying eqn. (3.4) by z m and summing over 0 to oo, we get 
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00 


I* 


00 00 CO 

Z p m+ i (n) z m - Z Pm (n) . z n = X S Pm (n-1) z m + p c Z p m (n+1) z" 
m=0 m=0 m=0 m=0 


00 


-(MkX)+ ZPm(n)Z n 

m=0 

00 00 00 00 

1/z Z Pm+1 (n) z m+1 - Z Pm (n) z n = X Z Pm (n-1) z m + |i c z Pm (n+1) z m 
m=0 m=0 m=0 m=0 


00 


- (p k X) + Z Pm (n) z m 
m=0 

1/z [P(z, m, (n) - po (n)] - P(z, m,n) = - X P(z, m, n-1) + p c P(z, m, n+1) -(k+\x c ) P(z, m, n) 
(l-z)/z P(z, m, (n) - po (n)] - 1/z p 0 (n) = - X P(z, m, n-1) + p 0 P(z, m, n+1) -(X+pc) P(z,m,n) 
s P (z, m, n) + (X + (i c ) P(z, m, n) +X. P(z, m, n-1) - p c P(z, m, n+1) = 1/z p 0 (n) 

(s+ X + p c ) p(z, m, n) - X p(z, m, n- 1 ) + p c p(z, m, n+ 1 ) = 1 /z p 0 (n) 

-X p(z, m, n-1) + (s+ 1 + p c ) p(z, m, n) + p c p(z, m, n+l) = l/zp 0 (n) (3.8) 

AP = [5y 5ii Siw] 7 

Where A is a real tri-diagonal matrix of order (N+1 )x (N+1) matrix. 

P is a column vector. 

8jj is the Kronecker delta defined as 


Sij 


y* if H 

0 ifj*i 
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. CHAPTER -FOUR 

A DISCRETE TIME TRANSIENT SOLUTION OF INTERDEPENDENT 
TANDEN QUEUEING MODEL WITH ARRIVAL IN BULK 

ABSTRACT 

In this Chapter we study a discrete time transient solution of an interdependent 
queueing model having three services in series. The customers arrive at the first node in 
groups of size k, the arrival of cutomers and number of service completion at each node are 
correlated, following correlated poisson process. Services are provided singly. For this 
queueing model, the average number of customers in the network, and the probability that the 
network is idle are derived. This model may be useful in finding the system performance 
measures for communication networks. 


INDRODUCTION : 

The communication system can be modelled as network of interconnected queues by taking 
the messages as customers communication buffer as waiting line and all activities necessary 
in transmission of the messages as services. In communication networks it is generally 
known that the packet switching gives better utilization over the circuit or massage 
switching and yields relatively short network delay. In packet switching the messages are 
divided in to small packets of random length each packet will have an independent header 
for routing. Due to the unproductive nature of demand, transmission lines congestion occurs 
in communication systems. This leads to model communication networks as a tandem 
queueing system. Several authors have studied the communication networks. Yakuo (1993) 
did throughout analysis of tandem queues, Jenq (1994) provided approximation for 


packetized voice traffic in statistical multiplexer in process but most of them have 
considered the independent assumption on the service and arrival process. However, in 
some communication systems like, store and forward communication network standard 
type of independence assumption is realistically inappropriate due to the fact that message 
generally preserve their length as they transverse the network. The inter-arrival and service 
requirements at queues internal to the system are thus dependent as there are queueing 
processes at each of the nodes of the network through which the same messages are routed. 
These dependencies can have a marked effect on system performance and must be 
accounted for any realistic analysis.Recently, Srinivasa et. al. (2000) worked on 
interdependent tandem queueing model, who developed and analyzed an interdependent 
communication network which is composed of several nodes of transmission through direct 
dependents between arrivals and transmission. But, still there seems to be no work done on 
the interdependent tandem queueing models with bulk arrive, although with reference to 
communication networks, it is most natural with respect to actual operations of such systems. 

It is in the reference to the communication networks, is assume that the massages are 
packetized at source stored in buffer for transmission. So for the concerned queueing model 
we have considered that customers are arriving in batches of size K. There are three nodes 
(servers) in series, after being served at the first node the customer leaves the network with 
probability pf or join the second buffer with probability pi after being transmitted to the 
second node the customer join the third buffer with probability (1 - pi) and leaves the system 
with probability (1 - P 2 ). 
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Here we assume that the arrival of customers (Packetes) and number of transmission 
services at each node (server) are correlated and follow correlated poisson process having the 
joint probability mass functions of the form (parallel to Rao etal. 2000). 

Here we assume that the arrival of customers (Packetes) and number of transmission 
services at each node (server) are correlated and follow correlated poisson process having the 
joint probability mass functions of the form (parallel to Rao etal. 2000). 

min (x b x 2 , x 3 , X 4 ) 

P (X, =x,, X 2 , = x 2 , X 3 = x 3 , X 4 = x 4 ; m) = e*^ 1 + m2 + m3 + - 3e >m I 

1=0 

[kX -e)m ]Y j [pi -e)m ]Y [p.2 -e)m ]Y J [p 2 -e)m ]Y J (€my 

(xi -j)! (x 2 -j)! (x 3 -j)! 

Xi,X 2 ,X 3 ,X4=0, 1,2.............. 

0< (1, p, s) < min (X, pi), i = 1, 2, 3 

Here 

kX = mean arrival rate of packets 

pii = mean transmission rate of first node 

(i 2 = mean transmission rate of second node 

p 3 = mean transmission rate of third node 

6 = mean dependence rate of the communication net work. 

The schematic diagram of network is given, 

Steady state difference differential Equations of Network : 

Let P m (m, n 2 , n 3 ) be the probability that there are ni packets in the first buffer, n 2 

packets at the second buffer, n 3 packets at the third buffer at time t. Then the difference - 

differential equations of the communications network with three node are : 

P m+ i (m, n 2 , n 3 ) - P m (m, n 2 , n 3 ) = - (X + pi + p 2 + p 3 - 4e) P m (ni, n 2 , n 3 ) + (kX-e) 

P m (n r k, n 2 , n 3 ) + (pi - €) pi P m (m+1, n 2 , n 3 ) + Pm(m+1, n 2 -l, n 3 ) + (p 2 - e) pi 

105 ; 















Pm(ni+1, n 2 , n 3 ) + ((J-i - e) pi (1 -p 2 ) P m (ni, n 2 +l, n 3 -l) + (p 2 - e) p 2 + P m (ni+1, n 2 , n 3 +l) ' j 

for ni, n 2 , n 3 > 0 .... (4.1) 

Pm+i (0, n 2 , n 3 ) ~ Pm (0, n 2 , n 3 ) - - (A + p 2 + p 3 - 3 e) P m (0, n 2 , n 3 ) + (pi - e) pi j 

Pm(l , n 2 , n 3 ) + (p 2 - e) (1 - p 2 ) P m (ni+1 , n 2 - 1, n 3 ) + (p 2 - e) p 2 P m (0, n 2 +l , n 3 ) + 

jfl 

(P 2 . - e) (1 - p 2 ) P m (l, n 2 +l, n 3 -l) (p 3 - e) P m (0, n 2 , n 3 +l) + for 0, n 2 , n 3 > 0. (4.2) |j 

Pm+i(ni, 0,n 3 )-P m (ni, 0,n 3 )--(A, + pi + ja 3 -3e)P m (ni, 0,n 3 ) + (kX-e)P m (m-k, 0,n 3 )+ ' 

(pi - e) Pi Pm(m+1, 0, n 3 ) + (p 2 - e) p 2 P m (m, 1, n 3 ) + (p 2 - e) (1 - p 2 ) P m (ni, 1, n 3 -l) + (p 2 - 

e) (1 -p 3 )P m (l, 0, n 3 -1) for m, n 2 > 0, n 2 = 0 ,...(4.3) ! 

P m +i (ni, n 2 , 0)-P m (ni, n 2 , 0) = - (A + pi + p 2 -3 e) P m (ni,n 2 ,0) + (kA -e) P m (nj-k, n 2 , 0) 

(pi - e) pi P m (ni+1, n 2 , 0) + (pi - e)(l - pi) P m (ni+l,n 2j 0) + (|x 2 - e)p 2 P m (ni, n 2 +l,0)+ | 

(p 3 - e)(l -p 2 )P m (ni, n 2 , 1) for m, n 2 , > 0, n 3 = 0 ...............(4.4) 

Pm+i (0, 0, n 3 )-P ra (0, 0, n 3 ) = - (A + p 3 - 2 g ) P m (0, 0, n 3 ) + (pi - e) pi P m (l, 0, n 3 ) + (jp 2 - e) 
p 2 P m (0, 1, n 3 ) + (p 2 - e) (1 -p 2 )P m (0, 1, n 3 -l) + (p 3 - e) P m (0, 0, n 3 +l) 

• for ni = n 2 = 0, n 3 > 0 ...............(4.5) ; ! 

P m +i (ni, 0, 0) — Pm (ni, 0, 0) = - (A + pi -2 g) P m (ni, 0, 0) + (kA -e) P m (ni-k, 0, 0) . 

(pi - e) pi P m (ni+1, 0, 0) + (p 2 - e) p 2 P m (ni, 1, 0) + (p 2 - e) (1 - p 2 )P m (ni, 1, 0) + 

(p 3 - e) P m (ni, 0, 1) for n 2 , n 3 , = 0, ni > 0 ........... . . ,.(4.6) || 

Pm+i (0, n 2 , 0) - Pm (0, n 2 , 0) = - (A + p 2 -2e) P m (0, n 2 , 0) + (pi - e) pi P m (l, n 2 , 0) + 

(pi - e) (1 -p 2 ) P m (0, n 2 -1, 0) + (p 2 - e) p 2 P m (0, n 2 + 1, 0) + (p 2 - e) P m (0, n 2 , 1) 

forni =n 3 , = 0,n 2 >0 ........(4.7) ■ fi 

P m +i (0, 0, 0) - Pm (0, 0, 0) = - (A - e) Pm(0, 0, 0) + (Pi - e) P m (l, 0, 0) + (p 2 - e) p 2 : 

P m (0, 1,0) + (p 3 - e)P m (0,0, 1) forni = n 2 = n 3 >0 ........(4.8) 

' . p§jj 

; v, 

' ' ' . ■ • . . ■ . ' ' ’ . .. . . ■ ■ ' . ■ iJl 
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Assuming the system is in steady state and using the boundary conditions and solving 
iteratively, we obtain the joint probability mass function of the contents in three buffers after 
reaching steady state as : 

a-e)V n 2 +n 3(l-p)V n 3(l-p) n 3 

Pm(ni, n 2 , n 3 ) = 

(pi - e) n i (p. 2 - e) n 2 (p 3 - e) n 3 

(Hi - kX.) [(m - pie -kX (1 - p)j [p 3 + { Pl , p 2 - (pi+ p 2 }.s - k k (1- pQ (1 - p 2 >] 

(pi - s) (^2 - e) (p 3 - e) 

(p 2 -s) (H2-S) 

for (k - s) < minimum[(pi - e), , ] 

(1 -Pi) (1 -Pi) (4-9) 

From this equation we can analyze the system performance measures by obtaining 
the system characteristics. 

The probability that the network is idle is : 

P m(0, 0, 0) = 

(Pi - kX) [(p 2 - Pl€ -kX (1 - pQ] [p 3 + { Pi- P2 - ( P1+ P2}-s - kX (1- Pi) (1 - p2 )] 

(pi - E) (|I 2 - E) (p 3 - s) 

(4.10) 

From this equation we see that probability of the idleness of the network increases 
when the mean dependent rate increases for fixed value of the other parameters : 

The joint probability of having ni packets in the first buffer and n 2 packets in the second 
buffer is 

(X-.€)V n 2 (l-p) n 2 

Pm(hl, ft2 5 ) — 

(^-e) n i(p2-e) n 2 

(pi - kk) [(pi - pie -kk (1 - pi)] 

(pi-8)(p 2 -£) (4.11) 


The joint probability of having n 2 packets in the second buffer and ns packets in the 
third buffer is : 


Pm(n 2 , n 3 ) 


(l-e)V n 3(l-pi)V n 3(l -p) n 3 

(P2 - e) n 2 (H3 - e) n 3 


(pi -kX) [(H2-Pie -kfr(T- Pi)] [P3 + { Pi. p 2 -(pi+P2}-g-kX(l-pi)(l p 2 )] 

(pi - 8) (p 2 - 8) (p 3 - e) 

( 4 . 12 ) 

The joint probability of having xi packets in the second buffer and x 3 packets in the 
third buffer is : 


Pm(ni, n 2 , n 3 ) = 


(X. - e) n , +n 3 (1 - pi) n 3 (1 - 
(Hi - e) n , (p 3 - e) n 3 

0*1 - kX) [(p 3 +{ Pi. P2 -(Pl+P2}.S-kX (1— pi) (1 — p2>] 


Oil - £ ) 0*3 - s) 


.( 4 . 13 ) 


The average number of packets in the network is : 

L = , ' ■ 

(A-i -s) [{p 2 -pie -kX(l — Pi) {p 3 + s{ pi.p 2 -(pi+ p 2 } -kX (l-pi)(l - p 2 )} 

+ (1 - pi) Oil - kX) {H3 + s{ Pi- P 2 - ( pi+ P 2 } - kX (1- pi) (1 - p 2 )} 

+ (1 - pQ (1 - p 2 ) Qii - kX) {p2 - kA. (1- pQ}] . 

Oil - kX) {|i 2 -Pie -kX (1 - p) {p 3 + s{ pi. p 2 - ( pi+ P 2 > — kX (1- pi) (1 - P 2 )} 

From this equation we observe that the mean number of packets in the network is a an 
increasing function of X. The average buffer content of the network is a decreasing function 
of mean transmission rates at first, second and third nodes when other parameters are fixed. 
The average buffer content of the network is a decreasing function of the mean dependence 
rate for fixed values of the other parameters. So by regulating the dependence parameter, one 
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can increase or decrease the mean number of packets in all the buffers and in each : 

intermediate buffer. The joining probabilities of the network 1 - pi and 1 - p 2 have an 
influence on the mean number of packets in the buffers. By decreasing 1 - p 1 and 1 - p 2 one : 

can reduce the conquestion in the intermediate queues. 

The mean number of packets in first buffer is : 

(X-6) 2 

Lqi = 

(m-e)(pi-€) 

Then mean number of packets in second buffer is : 

(>.-r) 2 n-p,/ 

Lqi = ■ 

(m - e) (pi - Pi s - kX(l - pi)) 

Then mean number of packets in third buffer is : 

<> - > >- <! -p.r 

Lq 2 = 

(M-3 - e) [p 3 + e{pi. P2 - ( P1+ P2>} - kx (1- pi)] 

The through put of the first transmitter is : 

m(^-e) 

m(l.-Po )= — 

pi (pi - kX) 

The average delay in the first buffer is : 

Li (m - e) 

pi(l-Po) pi(pi-kX) 

Through put of the second transmitter is : 

p 2 (^-e)(l-pi) 

1 * 2(1 -Po ) = 

* (H2 - e) 


109 




The average delay in the second buffer is : 

L2 (H2 - e) 

P 2 (l - Po) {2.2 [{22 - PlS - kX, (1 - Pi)] 


The through put of the third transmitter is : 

|i 3 (X-e)(l-pi)(l-p2) 

{2 3 (1 — P 0 .- •'*••) = 

({23 - e ) 

The average delay in the third buffer is : 

L 3 = (P3 - e) 

{l 3 (l -P 0 ...) {23 [{23 + e{pi.p 2 - (P1+P2)} - kX (1- Pi)] 
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A DISCRETE TIME TRANSIENT ANALYSIS OF A PRIORITY 
TANDEM QUEUE WITH ARRIVAL IN BATCHES AND HAVING NO 

INTERMEDIATE BUFFER . ; : . ■ 

ABSTRACT ' j,1 

In this chapter we have studied discrete time transient analysis priority tandem j 

queueing system with two stations in series. K type of customers arrive in batches of 
random size when all the customers in the batch belong to the same priority class and 
interarrival time of batches are distributed as exponentially. There is no intermediate buffer 

' ' , - | 

between the two stations, which may cause blocking at the first station. 

. ;■ . . > , 

Customers are served in a non-preemptive priority order. The expected delay in the 
system for the first member of a group of each type of customers is obtained, where all the 
customers have the same service time distribution. These type of queueing models are 
important in modelling and analysis of manufacturing systems and computer and 
co;onnuncation networks. 

INTRODUCTION: j 

. • ' . • ' ' ' ' ' ' ' ' ■' . • ''.'H 

. • ' . • . ■ ■ ■ . : / . ■ , ! 4 

Priority queues have important applications in modelling and analysis of m 

■ ■ ' .• • ' ’ ■ ■ . ;■ ' • . ■ 

manufacturing systems, computer systems and communication networks. Although, a lot 

of work has been done on priority queues, but, less attention is paid on priority tandem fg 

' ’ . . • ■ ' ; ■ •.*' , " ■. . * 

• . : ■ . . . . . m 

queues. : : ■ 9 
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In this chapter we consider a discrete time transient analysis tandem priority 

. 

' . ' ■" ’ .. . . ; ' ' ' . i| 

queueing system with two stations in series having no buffer between the stations, and K 

type of customers arriving in batches at the system according to discrete time 
homogeneous poisson processes. The lack of buffer between stations may cause blocking, 
because a customer finishing its service at the first station cannot move to the second 1 

station, if it is busy. Therefore the first station cannot start serving other customers until 1 

there is service completion at the second station. 


There has been an increased interest in tandem queues with blocking because of 
their wide applicability, see for example, Dattatreya (1978), Jun & Perros (1990), 
Xinltchao (1994) Grassman & Dekic (2000), However, no work has been done on priority 
tandom queue with blocking, having arrival in bulk. The model considered in this chapter 
is very common in communication system. For example in telecommunication system, 
message, are transmitted through the network. There is a finite buffer at each station and 
two types of message regular messages and acknowledgement messages. 
Acknowledgements messages carry information, about the stations of the system and 
have priority over regular messages. In packet switching network the message are 


decomposed into smaller pices called packet each of which has a maximum length so a 
message can be throught as a group of customers, when packets represent the customers. 


There is extensive literature on priority queues see Jaiswal (1968), Kleinrock 
(1976) and Wolff(1989) for extensive treatment and reference. 
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MODEL DESCRIPTION: 

In this model, we consider two station in series, with no intermediate buffer. We 
assume that customers arrive in batches of random size according to a discrete time 
homogeneous poisson process at the first station, with arrival rate (Number of 
groups/unit) for priority i customers (i = 1, 2 .... k). It is supposed that all customers in a 
batch classes are independent of each other. At station one, customers are surved singly in a 
non preemptive priority order with type i batch having priority over type j batch if i < j, 
within the same type the order of services is FCFS, within the same group. The order of 
services in random order i.e. batches of same priority are served in order of their arrival. 

When a customer’s service is completed station one it moves to station two 
provided station two is empty at that time i.e. arrival to station two is in single, if station 
two is busy the customer will wait at station one and block it until station two 
empty no customer can enter service at station one when it is being blocked. 

We assume that the service requirement of type i customers at the first station are 
i.i.d., r.v’s Si with distribution function F m (i) and density function f m (i). The station two all 
the customers have service time requirement T which are i.i.d., r.v's with distribution 
function G m . and density function g m - 

From the system structure, it is obvious that the order in which the customers are 
served at the first station is the same as that of the second station. The amount of time a 
customer spends at station two is exactly its service time. Hence if we can compute the 
expected time a customer spends at station one, we know the expected amount of time the 
customer spends in the whole system. 


The batch size X; is a random variable representing the number of customer in the 
group of priority i customer. 

P(xi = n) = gi, n i = l,2,......k, n= 1,2.. 

E(X,) = &\ and a is mean number of customers in an arbitrary batch i.e. a = E(a[). 

k 

X = Y Xi and 
i = 1 

00 

L S §i,n 

n= 1 

we define Pi = 

X 

ANALYSIS : 

We define the first station as busy if either, the server is serving a customer or is 
being blocked. If a customer I arrives at a busy station, the time from initiation of its 
service at station one to the moment of its departure to station which also called the 
occupation time of a customer at station one is : 

Uj = max {sj , T} 

because all customers have i.i.d. service requirements at the second station the 
■V. occupation times of first customer at the first station are independent. How ever, type i 

customer arrives at an empty station one the amount of time the customer stays their 
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dependents on the remaining service time at station two and its own service requirements 
at the first station. 


Therefore if only the first station is concerned it is exactly M x /G/1 priority 
queueing system with exceptional service time for the first customers of each busy period. 


Now we employ the work conservation approach to analyse the average delays in 
the system. We have assume that S(i) has distribution function F m (i) and density function 
f m (i) but it a batch of type i arrives an empty system, service requirements of its first 
customer is different from that of a regular customers in has the same distribution as Tjo 

with distribution function F m (iO) and density function f m (iO) where i = 1, 2 ,k. Therefore 

the system under consideration is be having like in M x /G/1 queue at station one with a 
compound poisson arrival process of rate X an in exceptional service time for customers 
who arrive at an empty system, because in order to find the average amount of work in a 
system we treat all the customers as being at the same type (With probability Pi = X-JX) 
which is a type i customers, now consider an arbitrary batch who arrives at a busy system 
the service requirement T of the first customer of this batch has distribution function F(m) 
and density function f(m) where 

k ?wSi 

F m = . £ F m (i) Pi = 

i = 1 Xa 

k 

F m = Z F m (i) 

i=l 


k k 

hence E(T) = Z Pi E(T io ), E(T 2 ) = ZPiE(T io 2 ), (5.1) 

i = 1 i=l 


Is a customer arrives at an empty system, he receives an exceptional service 
time To with distribution function F m (0) and density function f m (0), where 


k k 

F m (0) I Pi F m (io), F m (0) = Z Pi f m (0) 

i = 1 i=T 

k k 

Hence E(T 0 ) = Z Pi E(T io ), E(T 0 2 ) = Z Pi E(T io 2 ) ........(5.2) 

i= 1 i=l 

The average amount of work in the system, which is also average delay for an 
arbitrary customer is given by : 

XaE(T 0 2 ) + [k 2 aE(T 0 )E(T 2 )/ 1- *.aE(T)] 

E(V) = • • ....(5.3) 

2[l-XaE(T) + A,aE(T 0 )] 

Where E(T), E(T 2 ), E(T 0 ) and E(T 0 2 ) are given by (1) and (2). The proportion of 
time the server is idle then Po is given by : 

l-XaE(T) 

P 0 = - ..-(5.4) 

1 +A,aE(T 0 ) - XaE(T) 


Now we assume that Si is the service time of an arbitrary type i customer and S is 
the service time of an arbitrary customer, then 

E(S) = Z Pi E(S0, E(Si) = Po E(T i0 ) + (1 - P 0 ) E (Tj) 

E(Si 2 ) = Po E(T i0 2 ) + (1 - Po) E (T; 2 ), I = 1 , 2, ........ .k, 

The average amount of work in service is : 


k 

I h gi Po E(T i0 2 ) + (1 - P 0 ) E (Tj 2 )/2 
i= 1 

Where Po is given in eq.n (4) ( 5 . 5 ) 

We first analyse, Type 1 customers, the delay of the first customer of type 1 batch is 
the sum of the service times of types of type 1 customers waiting in the queue, and the 
remaining service time of any customer a type 1 batch, then 
di=Xia,d 1 E[T,] + XaE[S 2 ]/2 
= pidi + LaE [s 2 ] / 2, 

Where pi = Liaidi [Ti] therefore 

di = eidi + e 2 d 2 + A.a E [s 2 ] /2 

Then we analyse Type 2 customers. Let R be the amount of work in the system 
found by a randomly selected customer from an arriving Type 2 batch, the R is the sum of 
service times of all Types 1 and 2 batches waiting in the queue, plus the amount of work 
currently in service. Elence 

E (R) = LiaA E [T i] + L 2 a 2 d 2 E [T 2 ] + La E [s 2 ]/2 
= pidi + p 2 d 2 + La E [s 2 ]/2 

The expected delay of first customer of type 2 batches is the sum of E [R] and the 
expected service times of new Type 1 arrivals before the beginning of their service, which 
is d 2 LiaiE[Ti], 

d2 = E[R] + d2LiaiE[Ti] = E[R] + pidi 

pidi + LaE [S 2 ]/2 LaE(S 2 ) 

hence d 2 = — — = 

1-pi-pi 2(l-pi) (I-P 1 -P 2 ) 

■ ' W ' ■ 



I Mi[p 0 E(T io V(l-Po)E(Ti 2 ) 

i = 1 

2(1 - pi) (1 - pi - p 2 ) 

In general, let di be the average delay of the customer of a batch of i then, 


we have 


1/2 E(S 2 ) + E h } aj E[Tj] 
i = 1 

hence dj = 1 — - , i — 1,2 

1 -Gi 


Where G; - Z Pi 

j =1 

Therefore, the average delay for a type i customer is obtained as : 


XaE [S 2 ] 

2(1 - Gj.i) (1-Gj) 


I Wi ai [Po E(Tio 2 ) + (1 - P 0 ) E (Tj 2 ) 
i=l 

2(1 -pi) (1 - P 1 -P 2 ) 


Where pi = h a; E(Ti), and P 0 is given in (4). 

We define the first station as busy if other the server is serving a customer, or a 
customer is being blocked. If a batch of i customer arrive at a busy system, their for the first 
customer of the batch. The time from the initiation of its service at station one to the 
initiation of its service at station two, will be called as the occupation time of this customer 


at station one, which 


Ui - max {si , T} ; customers 



Fm - fm(i) G(m) + F m (i) g(m), i = 1, 2, k. 

However, if a group of i customer arrives at an empty station one, the amount of 
time the first customer of the batch stays there depends on the remaining service time at 
station two (if any) and its own service requirement at the first station. 

Assume that the first customer of the arrives at the first station units of time after 
the previous customer departs from the first station, and is a type i customer, then, the 
amount of time this customer stays at station one in max (Si, T - s}. Hence the conditional 
density function of the amount of time the first customer of type i stays in station one is 

F m (i) G(m+s) + f m (i) g(m+s). 

The probability density that a customer arrives at time s and is of class i is 

k 

X\ e x \ where again X = Y.X\. Therefore, given that a type i customer arrives at ah 

i = 1 

empty station one, the density function of the amount of time this customer stays in station 


one is 


F m (io) = f m (i) G(m) + e xl (X [F m (i) + f m (i)] Je“ Xs g(s) ds 


.(5.8) 


Now, We have obtained that the expected delay of the first customer of a type I 
batch at the first station is dj 

XaE [S 2 ] 

2(1-1 Pi) (1 -Z Pi) 


Where pi = X- t E[Ui], XE[S 2 ] = Z ^ a ; [P 0 E(U io 2 ) + (1 - p) E(Ui 2 ) and 

i=i . : ■ 

E[Ui], E[U 2 ], E[Ui 0 ] and E[Ujo 2 ] are computed from (5.7) and Po is given by 


1 - 1 Xi a; E[Ui] 

i=l 


k k 

1 + Z Xi ai E[U i0 ] + S Xi a; E[U;] 
i=l i=l 

The expected time, the first customer of a type i batch spends at the first station is 
di + P 0 E[U i0 ] + (1-Po) E[Ui], i = 

and the expected total time, the first customer , of a type i batch spends in the whole system 


d’i = di + P 0 E[U i0 ] + (1-Po) E[Uj] + E[T], i = 1, k, 


.(5.10) 


Where E[T] is the common mean of the service requirement at station two. From little’s 
formula, the expected number of type i customer in the system 

Li = A.j a, dj = Xi aj [d; + PoE(Ui 0 ) + (1-Po) E(Uj) + E[T], i = 1, ........k 
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